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Abstract 

Exotic atoms provide a unique laboratory for studying strong interactions and nuclear medium 
effects at zero kinetic energy. Experimental and theoretical developments of the last decade in 
the study of exotic atoms and some related low-energy reactions are reviewed. The exotic atoms 
considered are of vr", K~ , p, S~, and also the so far unobserved H~ atoms. The analysis of these 
atomic systems consists of fitting density dependent optical potentials Fopt = tip)p to compre- 
hensive sets of data of strong-interaction level shifts, widths and yields across the periodic table. 
These provide information on the in- medium hadron-nucleon t matrix t{p) over a wide range of 
densities up to central nuclear densities. For pions the review focuses on the extraction of the ttN 
in-medium s-wave interaction from pionic atoms, which include also the recently observed at GSI 
deeply bound vr" atomic states in isotopes of Sn and Pb. Also included are recent measurements 
at PSI of elastic scattering of vr^ on Si, Ca, Ni and Zr at 21.5 MeV. The experimental results 
are analyzed in the context of chirally motivated vr-nuclear potentials, and the evidence for partial 
restoration of chiral symmetry in dense nuclear matter is critically discussed. For antikaons we 
review the evidence from K~ atoms, and also from low-energy K~p scattering and reaction data 
for and against a deep ^-nucleus potential of 150-200 MeV attraction at nuclear matter density. 
The case for relatively narrow deeply bound K~ atomic states is made, essentially independent 
of the potential-depth issue. Recent experimental suggestions from KEK and DA<I>NE (Frascati) 
for signals of ^-nuclear deeply bound states are reviewed, and dynamical models for calculating 
binding energies and widths of K-nuclear states are discussed. For kaons we review the evidence, 
from total and reaction cross section measurements at the AGS (BNL) on Li, C, Si and Ca 
at piab = 500 — 700 MeV/c, for significant absorptivity of tKN{p) beyond that expected from ijpjy 
within the impulse approximation. Attempts to explain the extra absorptivity for the relatively 
weak interaction of K mesons in terms of a hypothetical exotic S = +\ pentaquark strength 
are reviewed. For antiprotons the exceptionally broad data base due to the recent results of the 
PS209 collaboration at CERN are analyzed, together with results of radiochemical experiments. 
We discuss the dependence of the phenomenological |)-nucleus interaction on the model adopted for 
the neutron density, showing how the neutron densities favored by our comprehensive analysis are 
compatible with densities from other sources, including our own analysis of pionic atoms. It is also 
shown how the strong absorptivity of the p- nucleus interaction, which leads to the prediction of 
saturation of widths in deeply-bound p- atom states, also explains the observed saturation effects 
in low-energy p annihilation on nuclei. For S hyperons we review the evidence, from continuum 
hypernuclear (vr", K^) spectra obtained recently at KEK on C, Si, Ni, In and Bi, for substan- 
tial repulsion in the S-nucleus interaction, and the relationship to the inner repulsion established 
earlier from the density-dependence analysis of S~ atoms and by analyses of past {K~ ,t:^) AGS 
experiments. 
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Lastly, for H hyperons we review prospects of measuring X-ray spectra in H~ atoms and thereby 
extracting meaningful information on the H-nucleus interaction. The significance of the latter 
to the physics of AA hypernuclei and to extrapolation into multistrange hypernuclei are briefly 
reviewed. 
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I. INTRODUCTION 
A. Preview 

In 1997 we published together with Chris Batty a Physics Reports review Strong Inter- 
action Physics from Hadronic Atoms [l[ that has had a substantial impact on the progress 
made in the last decade on the study and understanding of in-medium nuclear interactions 
for various hadrons at low energy. The most spectacular advance on the experimental side 
has been perhaps the discovery and study of deeply bound pionic-atom states using the 
'recoil-free' (d, '^He) reaction near the pion production threshold on isotopes of Sn and Pb 
at GSI. Another potentially promising advance, very recently, concerns the as yet weak ev- 
idence for the existence of deeply bound antikaon-nuclear states gathered from stopped K~ 
reactions studied at KEK and at DA$NE, Frascati. These advances were triggered by, and 
have stimulated related theoretical work in which low-energy in-medium hadronic properties 
were considered in terms of a systematic chiral-perturbation expansion. This holds not only 
for pions, owing to their relatively small mass, but also for antikaons where the dominant 
effect of the 5* = — 1 subthreshold quasi-bound state A(1405) was treated by a unitarized 
coupled-channel approach based on chiral-perturbation expansion. Significant advances in 
the study of nuclear interactions of other hadrons at low energy, for kaons, antiprotons and 
S hyperons, have also been made. A good evidence for the wide interest within the nuclear 
physics and hadronic physics communities in these subjects is provided by the large number 
of recent topical reviews devoted to the overall theme of in-medium nuclear interactions of 
hadrons, including energy regimes higher than considered here. A representative list of such 
Reviews during the last five years covers the following subjects: chiral symmetry in nuclei 
and dense nuclear matter 0], pions in nuclei, aprobe of chiral symmetry restoration chi- 
ral symmetry and strangeness at SIS energies [4], medium modifications of hadrons - recent 
experimental results [5|, kaon production in heavy ion reactions at intermediate energies 6|, 
nucleon and hadron structure changes in the nuclear medium and impact on observables 7|. 

In the present Review we discuss and summarize the developments in understanding 
the in-medium properties of several hadron-nucleon systems at low energy as unraveled by 
our recent phenomenological studies and related ones by other authors in this field. This 
brief Preview subsection is followed by brief introductory subsections on wave equations and 
optical potentials, on nuclear densities, and on in-medium interactions. 
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B. Wave equations and optical potentials 



The interaction of hadrons in nuclear medium of density p is traditionally described by 
a dispersion relation based on the Klein-Gordon (KG) equation 

E''-Y>^-m'-Il[E,Y>.p) = Q. n = 2£;Kpt, (1) 

where n(£', p, p) is the hadron self-energy, or polarization operator and Vopt is the optical 
potential of the hadron in the medium [Sj . Here m, p and E are the rest mass of the hadron, 
its three-momentum and energy, respectively. For finite nuclei, and at or near threshold as 
applicable to most exotic-atom applications, Eq. ([1]) gives rise to the following KG equation: 

V2 - 2p{B + Kpt + K) + (K + Bf] V; = (/i = c = 1) (2) 

where /i is the hadron-nucleus reduced mass, B is the complex binding energy and 
is the finite-size Coulomb interaction of the hadron with the nucleus, including vacuum- 
polarization terms, added according to the minimal substitution principle E ^ E — Vc- A 
term 2V^V^pt and a term 2BVopt were neglected in Eq. ([2]) with respect to 2//V^pt; the term 
2SV^pt has to be reinstated in studies of deeply-bound states. 

The simplest class of optical potentials Vopt is the generic tp{r) potential, which for 
underlying s-wave hadron-nucleon interactions assumes the form: 

2/iKpt(r) = -47r(l + {bo[pn{r) + Pp{r)] + rMpnir) - Pp{r)]} . (3) 

Here, p„ and Pp are the neutron and proton density distributions normalized to the number 
of neutrons and number of protons Z, respectively, M is the mass of the nucleon and 
Tz = +1 for the negatively charged hadrons considered in the present Review^. In the 
impulse approximation, bo and hi are minus the hadron-nucleon isoscalar and isovector 
scattering lengths, respectively. Generally these 'one-nucleon' parameters are functions of 
the density p, but often the density dependence may be approximated by fitting effective 
values for bo and bi to low-energy data. Extensions to situations which require 'two-nucleon' 
terms representing absorption and dispersion on pairs of nucleons, or which are motivated 
by p-wave hadron-nucleon interactions, will be dealt with in the next section, for pionic 
atoms. 

For scattering problems, the applicable form of the KG equation is given by: 

V' + e- {2eitKv, + Kpt) - K')] ^ = (4) 

in units of h = c = 1, where k is the wave number in the center-of-mass (cm.) system. For 
the simplest possible tp s-wave term, the optical potential l^pt is of the form 

242Kpt(r) = -4nFA{bo[pn{r) + Pp(r)] + rMPn{r) - Pp(r)]} , (5) 



where e^f^ is the cm. reduced energy. 



[e 



red / -'-^p 



E;' + (6) 



Tz = —2tz for isodoublets and —t^ for isotriplets, where tz is the value of the z th projection of isospin 
for the hadron considered. 
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in terms of the cm. total energies Ep for the projectile and Ea for the target, and 

^ M{Ea + Ep) 



(7) 



is a kinematical factor resulting from the transformation of amplitudes between the hadron- 
nucleon and the hadron-nucleus cm. systems. Here Ma is the mass of the target nucleus and 
^/s the total projectile-nucleon energy in their cm. system. These forms of the potential and 
the equation take into account 1/A corrections, which is an important issue when handling 
light nuclear targets. The kinematical factor Fa reduces at threshold to the kinematical 
term (1 + (1 — 1/A)fi/M) appearing in Eq. Q for the one-nucleon s-wave potential term in 
exotic-atom applications. 

For Fermions, such as antiprotons or E hyperons, one might ask why the KG equation 
is here used in going into the relativistic domain instead of the Dirac equation. Indeed 
when interpreting experimental transition energies, in order to extract the strong interaction 
effects, it is essential to use the Dirac equation with finite size nuclear charge distribution and 
vacuum polarization terms, e.g. Ref. |9| for antiprotonic atoms. However, strong interaction 
effects are normally given as proper averages over the fine structure components. The use of 
the KG equation rather than the Dirac equation is numerically justified when fine-structure 
effects are negligible or are treated in an average way, as for the X-ray transitions considered 
here. The leading j dependence (j = / ± i) of the energy for solutions of the Dirac equation 
for a point-charge 1/r potential goes as (j + ^)~\ and on averaging it over the projections 
of j gives rise to (/ + |)~^ which is precisely the leading / dependence of the energy for 
solutions of the KG equation. The higher-order contributions to the spin-orbit splitting 
are suppressed by OlZa/n)"^ which, e.g., is of order 1% for the high-n X-ray transitions 
encountered for antiprotons. It was checked numerically for few typical cases that the spin- 
orbit averaged shifts and widths thus obtained differ by less than 1% from the (2j+l)-average 
of the corresponding quantities obtained by solving the Dirac equation. This difference is 
considerably smaller than the experimental errors placed on the measured X-ray transition 
energies and widths. 



C. Nuclear densities 

The nuclear densities are an essential ingredient of the optical potential. The density 
distribution of the protons is usually considered known as it is obtained from the nuclear 
charge distribution ^] by unfolding the finite size of the charge of the proton. The neutron 
distributions are, however, generally not known to sufficient accuracy. A host of different 
methods have been applied to the extraction of root-mean-square (rms) radii of neutron 



distributions in nuclei but the results are sometimes conflicting, e.g. Refs. [11|, [12, ll3|, 114 
11,0. For many nuclei there is no direct experimental information whatsoever on neutron 
densities and one must then rely on models. To complicate things further we note that there 
is a long history of conflict between values of neutron rms radii derived from experiments 
using hadronic projectiles and neutron rms radii obtained from theoretical calculations. For 
that reason we have adopted a semi-phenomenological approach that covers a broad range 
of possible neutron density distributions. 

Experience with pionic atoms showed [13] that the feature of neutron density distributions 
which is most relevant in determining strong interaction effects in pionic atoms is the radial 
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FIG. 1: Fit of a linear expression in the asymmetry parameter to RMF values of Vji — Vp. 



extent, as represented e.g. by r„, the neutron density rms radius. Other features such as 
the detailed shape of the distribution have only minor effect. For that reason we chose the 
rms radius as the prime parameter in the present study. Since r^, the rms radius for the 
proton density distribution, is considered to be known, we focus attention on values of the 
difference — rv A linear dependence of r„ — rp on (A^ — Z)/A has been employed in p 



studies [ly, [181, Il9[ , namely 

N-Z , 

rn-rp = 7 ^ +0 , (8) 

with 7 close to 1.0 fm and 6 close to zero. The same expression with 7 close to 1.5 fm was 
found [17[ to represent well the results of relativistic-mean- field (RMF) calculations Q for 
stable nuclei, as shown in Fig. [H but these values of r„ — are larger by about 0.05-0.10 fm 
than the 'experimental' values in medium-weight and heavy nuclei used in recent relativistic 
Hartree-Bogoliubov (RHB) versions of mean-field calculations 21, 2^- Expression ([8]) has 



been adopted in the present work and, for lack of better global information about neutron 
densities, the value of 7 was varied over a reasonable range in fitting to the data. This 
procedure is based on the expectation that for a large data set over the whole of the periodic 
table some local variations will cancel out and that an average behavior may be established. 
Phenomenological studies of in-medium nuclear interactions are based on such averages. 

In order to allow for possible differences in the shape of the neutron distribution, the 'skin' 
and 'halo' forms of Ref. [18] were used, as well as an average between the two. Adopting a 
two-parameter Fermi distribution both for the proton (unfolded from the charge distribution) 
and for the neutron density distributions 

(r) = , (9) 

"'^ 1 -h exp((r - Rn,p)/an,p) ' 
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FIG. 2: Proton and neutron finite-range folded densities for ^^''Sn with (3 = 0.85 fm. Neutron 
densities are calculated for 7 = 1.2 fm, see Eqs. ([8l [T0|) . 

then for each value of r„ — Vp in the 'skin' form the same diffuseness parameter for protons 
and neutrons, a„ = Op, was used and the Rn parameter was determined from the rms 
radius r„. In the 'halo' form the same radius parameter, Rn = Rp, was assumed and 
was determined from r„. In the 'average' option the diffuseness parameter was set to be 
the average of the above two diffuseness parameters, a^™ = (a^ + a^)/2, and the radius 
parameter Rn was then determined from the rms radius r„. In this way we have used three 
shapes of the neutron distribution for each value of its rms radius all along the periodic 
table. These shapes provide sufficient difference in order to be tested in global fits. 

Another sensitivity that may be checked in global fits is to the radial extension of the 
hadron-nucleon interaction when folded together with the nuclear density. The resultant 
'finite range' density is defined as 

p^ir) = /rfrV(r')-^e"(-^')V/^^ , (10) 

assuming a Gaussian interaction. Other forms such as a Yukawa function may also be used. 
Figure [2] shows for example finite-range folded proton and neutron densities in ^^''Sn, calcu- 
lated using the three models listed above for generating neutron densities. The difference 
between these three models becomes pronounced from about 8 fm on, a radial extent to 
which p atoms are particularly sensitive. 

D. In- medium interactions 

The tp form of the optical potential, where t is the two-body hadron-nucleon t matrix and 
p is the nuclear density (more precisely, the nucleon-center distribution density), holds at 
high collision energy where most of the nuclear medium effects such as the Pauli principle are 
negligible. At low energy, and particularly near threshold, nuclear medium effects may and 
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FIG. 3: Real (top) and imaginary (bottom) parts of the in-medium isospin-averaged K~ N (ef- 
fective) threshold scattering amplitude as function of density, calculated in Ref. [3] including the 
effects of the Pauli principle (dashed lines), plus the effect of self consistency for propagation 
(dot-dashed lines), and also the effect of self consistency for propagation (solid lines). 

often do assume special importance. To demonstrate the scope of medium effects we use as an 
example the case of K~ nuclear interaction near threshold {yfs ~ rnx- + Mp = 1432 MeV), 
for which the underlying two-body KN system is strongly coupled to the vrS and ttA reaction 
channels, all in s waves. In particular, the interaction in the 1 = KN-ttTi coupled-channel 
system is so strong as to generate a subthreshold quasibound KN state about 27 MeV below 
threshold with a width approximately 50 MeV. This A (1405) is observed as a resonance in 
ttS final-state interaction spectra ^]. The K~ nuclear optical potential in the large A limit, 
from Eq. ([3]), assumes the form: 

2/iKpt(r) = -47r(l + ■^)[aK'p{p)pp{r) + aK-n{p)pn{r)] , (11) 

where 

o-K-p = bo - bi , ax-n = bo + bi (12) 
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in the particle basis. The density dependent effective scattering amphtudes a^-pip), CLK-nip) 
are complex due to the coupling to the reaction channels. For p ^ 0, the low-density limit 
asserts their limiting values ax-p, clk-u respectively, where the latter are (strictly speaking 
minus) the corresponding scattering lengths. Figure [3] shows the density dependence of 
the effective isoscalar threshold scattering amplitude a^s = |[ax-p(p) + cbR-nip)] for three 
cases: (i) no medium effects beyond Pauli blocking are included (dashed line); (u) a self- 
consistent calculation including the K self energy (dot-dashed line); and (iii) including 
also the nucleon self energy (solid line). The change of the sign of Reaefr from negative to 
positive corresponds to the transition from an apparently repulsive free-space interaction to 
an attractive one in the nuclear medium. The underlying physics is that the Pauli principle 
suppresses the contribution from Pauli forbidden intermediate states, thus weakening the 
in-medium t matrix which no longer supports a subthreshold quasibound state; thus, the 
A(1405) gets pushed up above threshold 25, 2^. The inclusion of K and (self energy) 



medium modifications pushes this transition to a lower density, as first shown by Lutz [27|, 
but the free-space (p = 0) threshold scattering amplitude remains negative, reflecting the 
dominance of the A(1405) / = subthreshold resonance. The figure shows that apart from 
the very low density regime, the K~ optical potential Vopt = t{p)p evaluated within such 
self consistent models is well approximated over a wide range of densities by a tesp form, 
where tes = — (27r/yu)(l + p/M)aes = const. A genuine p dependence of tes appears only at 
very low densities. The strength of Reaeg is seen to be reduced to about 50% of its initial 
value upon imposing self consistency. This is due to the suppressive effect of Imtefr in the 
K~ propagator of the Lippmann-Schwinger equation for t{p): 

t = V + V 1. (13) 

E-HZ-tp-V^ + iO 

Here v and t are coupled-channel meson-baryon (mB) potential and t matrix, respectively, 
and H^^^ is the corresponding kinetic energy operator which depends implicitly on the density 
p through the imposition of the Pauli principle in KN intermediate states. The K~ optical 
potential tp and the nucleon potential Vn act only in KN intermediate states. A sizable 
Imt leads to an exponential decay of the propag ator {E - H^^l - tp - V/v + iO)~\ so that 
t ^ V thus losing the cooperative coupling effect to the irY channels in higher-order terms 
of V. 



II. EXPERIMENTAL BACKGROUND 

In this section we outline experimental results, relevant to the topic of the present work, 
that have been obtained since the publication of Ref. [I| which was focused on the strong 
interaction physics involved in exotic atoms of mostly medium weight and heavy nuclei. As 
before, most of the information on the interaction of low energy hadrons with nuclei which 
provides insight to in-medium properties, comes from strong interaction effects in exotic 
atoms. Here we include only a brief reminder of exotic atoms, referring to other reviews for 
more details [H, [2^. 

An exotic atom is formed when a negatively charged particle stops in a target and is 
captured by a target atom into an outer atomic orbit. It will then emit Auger electrons and 
characteristic X-rays whilst cascading down its own sequence of atomic levels until, at some 
state of low principal quantum number n, the particle is absorbed due to its interaction with 
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the nucleus. The hfetimes of all the particles considered here, namely vr", K~ , p and S~, are 
much longer than typical slowing down times and atomic time scales. Therefore, following 
the stopping of the hadron in matter, well-defined states of an exotic atom are established 
and the effects of the hadron-nucleus strong interaction can be studied. The overlap of the 
atomic orbitals with the nucleus covers a wide range of nuclear densities thus creating a 
unique source of information on the density dependence of the hadronic interaction. 

In the study of strong interaction effects in exotic atoms, the observables of interest are 
the shifts (e) and widths (F) of the atomic levels caused by the strong interaction with 
the nucleus. These levels are shifted and broadened relative to the electromagnetic case 
but the shifts and widths can usually only be measured directly for one, or possibly two 
levels in any particular exotic atom. The broadening due to the nuclear absorption usually 
terminates the atomic cascade at a low value of the radial quantum number n, thus limiting 
the experimentally observed X-ray spectrum. In some cases the width of the next higher n+1 
'upper' level can be obtained indirectly from measurements of the relative yields of X-rays 
when they depart from their purely electromagnetic values. Shifts and widths caused by the 
interaction with the nucleus may be calculated by adding an optical potential to the Coulomb 
interaction. The study of the strong interaction in exotic atoms thus becomes the study of 
this additional potential. On the experimental side, studies of strong interaction effects in 
exotic atoms have been transformed over the years with the introduction of increasingly 
more advanced X-ray detectors and with increasing the efficiency of stopping the hadrons, 
such as with a cyclotron trap [28.]. In recent years exotic atom physics has turned into 
precision science. 

With the present topic of in-medium interactions of low energy hadrons we include not 
only data on exotic atoms but also data on the interaction of hadrons with nuclei at low 
kinetic energies where the interaction models have similarities to the models used with 
exotic atoms. In such cases features of the interaction may be studied across threshold, 
thus enhancing our knowledge of the hadron-nucleus interaction. Moreover, the respective 
free hadron-nucleon interactions at very low energies are obviously the reference to which 
the in-medium interactions have to be compared. Therefore key experiments on those more 
elementary reference systems will also be mentioned in the present section. 

Starting with pions, recent years have seen the continuation of experiments at PSI on 
pionic hydrogen and pionic deuterium, with ever increasing sophistication and efficiency, 
reaching accuracies which are limited by theoretical corrections [1^. Pionic atoms of deu- 
terium are obviously the source of knowledge on the pion-neutron interaction, which is a 
pre-requisite for studies of heavier targets and for obtaining separately the isoscalar and the 
isovector interactions. The use of deuterium for this purpose inevitably introduces some 
dependence on models in the extraction of the two basic interactions. 

Turning to heavier pionic atoms, the last decade has been dominated by the experimental 
observation of 'deeply-bound' pionic atom states in the recoil-free (d,^He) reaction jsoj which 
populates such states from 'inside' the nucleus. With this technique one avoids the cut-off 
imposed by nuclear absorption on the usual process of X-ray emission during the atomic 
cascades. The observation of these states in isotopes of Pb and Sn [sH made it possible 
to test predictions of interaction models which were based on data for conventional pionic 
atom states. Alternatively, it became possible to derive interaction parameters from deeply 
bound atomic states, for comparisons with models based on X-ray data. Moreover, studies 
of several isotopes of a single element have the promise of providing information on the role 
played by neutron density distributions in the pion-nucleus interaction. 
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Motivated by the renewed interest in the pion-nucleus interaction at very low energy, 
caused by the observation of deeply-bound states and by the possibility of linking the long- 
standing s-wave 'anomaly' with aspects of chiral-symmetry restoration (see below), differ- 
ential cross sections for the elastic scattering of both vr^ and tt^ by several target nuclei 
were measured at 21 MeV. A dedicated experiment H, 33 1 where both charge states of the 
pion were measured with the same setup and where absolute normalizations were provided 
by muon scattering, yielded angular distributions which could be analyzed with the same 
interaction models as used for pionic atoms, thus providing tests across thresholds of various 
characteristics of the pion-nucleus interaction. 

Turning to kaonic atoms, significant progress has been made with kaonic hydrogen, due to 
experiments with precision greatly exceeding that of the earlier generation experiments that 
removed the so-called kaonic hydrogen 'puzzle', where the strong interaction shift appeared 
initially to be attractive, contrary to expectations. The experiments at KEK [34,] showed 
that the shift of the Is level in kaonic hydrogen is repulsive, as expected from earlier phase- 
shift analyses. More recent results from DA<I>NE [35|, using the unique low energy K~ and 
from decay, are barely consistent with the KEK results within error bars. In addition, 
experiments at KEK on kaonic atoms of ^He seem to produce results quite different 
from the previous ones, which were at variance with predictions of most calculations. With 
the exception of the above two examples, the world's data on kaonic atoms have not been 
expanded in the last decade, because low energy K~ beams of sufficient quality are not 
available. 



Outside the realm of kaonic atoms, there have been experimental indications |37l . l38l . [3; 



40[ | of possible existence of strongly bound antikaon states in nuclei. These caused renewed 



interest in the question of the depth of the real part of the /^"-nucleus potential at thresh- 
old, where 'deep' real potentials are known from fits to kaonic atom data and 'shallow' 
potentials are obtained from chiral approaches. Although the experimental situation is not 
settled at the present time, some understanding of the antikaon-nucleus interaction is being 
promoted thanks to studies inspired by the speculations on strongly bound states. 

The interaction of mesons with nuclei has not been discussed in our previous Re- 
view [l[. This topic is included in the present Review since further analyses of previous 



transmission experiments have clearly demonstrated [41|, |42|, |43| that the elementary KN 
interaction in the range of 500-700 MeV/c is modified in the nuclear medium. In common 
with the case of K~ mesons, the recent renewed interest in the topic of medium-modification 
of the elementary interaction was motivated by speculations based on experimental in- 
dications, in this case on the possible existence of the G"*" pentaquark. 

The large cross sections for annihilation of antiprotons on nucleons set the scene for 
the interaction of p with nuclei at low energies. In terms of optical potentials that means 
dominance of the imaginary part which complicates the issue of the connection between the 
free pN interaction and the interaction in the nuclear medium. It also means that p do not 
penetrate deeply into nuclei. On the other hand, one of the experimental consequences is 
the ability to measure annihilation cross sections at extremely low energies. Indeed such 
measurements have been made in recent years at momenta as low as 40-50 MeV/c on light 
nuclei 0, 45] and on hydrogen Total cross sections for np down to 50 MeV/c have 



also been reported |4j 

The experimental situation with p atoms has changed significantly in the last decade with 
the publication by the PS209 collaboration 4^ of high-quality data for several sequences of 
isotopes along the periodic table. For most target nuclei strong-interaction level shifts and 
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widths have been measured for two atomic levels of the same antiprotonic atom, with two 
examples where information is available for three levels. In addition close to 20 target nuclei 



have been studied by the radiochemical method [49|, [SOj, observing the production of nuclei 
differing from the target nucleus by the removal of one neutron or one proton. Such data 
provide unique information on the absorption of p by a neutron or by a proton, repectively, 
at about 2.5-3 fm outside of the nuclear surface, and in particular the ratios between the 
probabilities for the two processes are determined quite reliably. The two kinds of data, 
namely, level shifts and widths on the one hand and ratios from radiochemical data on the 
other, were shown to lead to consistent results and could also be analyzed together. 

Finally we mention the interaction of low energy S~ with nuclei. No additional data 
on Zl~ atoms have been produced in recent years and the only relevant new experimental 
information was obtained from the {tt~,K~^) reaction on nuclear targets |51,] which showed 
some features in common with the atom potential derived more than a decade earlier. 



III. PIONS 

A. The pion-nucleus potential 

At zero energy the interaction of pions with nucleons is rather weak and consequently 
a tp approach would be expected to yield a reasonably good optical potential or at least 
provide a theoretically motivated form for the potential. The interaction of low energy pions 
with nucleons is affected by the (3,3) resonance at about 180 MeV, and this is manifested 
by a significant p wave term in the vrA^ interaction which in turn is reflected in the form of 
the optical potential j8|. The Kisslinger potential, where this p wave interaction leads to 
gradient terms in the pion nucleus potential, was introduced more than half a cen tury ago 
5^ and model-independent analyses of elastic scattering of pions by nuclei showed jHsf that 



indeed the local-equivalent potential has all the features expected for the Kisslinger potential 



(see also jM]). A tp potential at zero energy is real because pions cannot be absorbed at 
rest by a single nucleon, although they can be absorbed by the nucleus. For that reason 
Ericson and Ericson j55| introduced p^ terms into the potential which describe schematically 
the absorption of tt" on pairs of nucleons. The potential for tt" mesons then becomes, in its 
simplest form, 

2/xKpt(r) = g(r) + V-a(r)V (14) 

with 

g(r) = -47r(l + ^){6o[p„(r) + pp(r)]+6i[p„(r)-pp(r)]} 

-M^ + ^HBoPn{r)pp{r), (15) 



a(r) = 47i{l + -^)-'{co[pn{r)+Pp{r)]+c,[pn{r)-pp{r)]} 

+4vr(l + ^)-^4Cop„(r)p,(r), (16) 

where p„ and Pp are the neutron and proton density distributions normalized to the number 
of neutrons N and number of protons Z, respectively, p is the pion-nucleus reduced mass 
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and M is the mass of the nucleoli. In this potential q{r) is referred to as the s wave potential 
term and a{r) is referred to as the p wave potential term. The real coefficients 60 and hi 
are minus the pion-nucleon isoscalar and isovector s wave scattering lengths, respectively, 
whilst the real coefficients Cq and Ci are the pion-nucleon isoscalar and isovector p wave 
scattering volumes, respectively. The parameters Bq and Cq represent s wave and p wave 
absorptions, respectively, and as such have imaginary parts. Dispersive real parts are found 
to play a role in pionic atom potentials. The PnPp in the absorption terms represent two- 
nucleon absorption which takes place predominantly on neutron-proton pairs. These were 
originally written jsHj as -BoPm Copfn, with pm = pn + Pp without distinguishing between 
neutrons and protons. The factor 4 is introduced above to make the coefficients Bq and Co 
comparable in the two formulations. In practice when parameters are obtained from fits to 
the data the two forms yield practically the same results. 

In the above expressions the terms linear in the nuclear densities are associated, in the tp 
approach, with the interaction between pions and free nucleons. Ericson and Ericson showed 
that the p wave dipole interaction is modified in the nuclear medium in a way analogous 
to the Lorentz-Lorenz effect in electrodynamics, replacing the above expression for a{r) as 
follows: 

where ^ is a constant of the order 1. This effect is generally referred to as the Lorentz- 
Lorenz-Ericson-Ericson (LLEE) effect and it results from short range repulsive correlations 
between nucleons. This modification should apply only to the linear part of a, and Eq. (|T6l) 
is re-written as 

with 

«i(r) = 47r(l + A)-i{co[p„(r) + p,(r)] + Ci[p„(r) - p,{r)]} (19) 

a2(r) = 47r(l + J^)-HCoPn{r)Pp{r). (20) 

Another complication arises due to the parameter 60 being exceptionally small. Hence 
second order effects in the construction of the isoscalar s wave potential term in q{r) become 
important jssj and this causes 60 to be replaced by 

bo = bo-^ibl + 2bj)kF, (21) 

where kp is the Fermi momentum calculated for the local nuclear density. 

Finally there is another relatively small term of a kinematical origin whose presence is 
supported by fits to pionic atom data. This is the so-called angle transformation term 
m |57j which is given by 



2/iAKpt = -47r{^(l + ^)-V[co(p„ + Pp) + ci(p„-pp)] 



+^(1 + ^ 



(22) 
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TABLE I: Data for Is states in pionic atoms 





shift (keV) 


width (keV) 


Ref. 


20Ne 


-32.17 ± 0.77 


15.43 ± 0.41 


[60] 


22Ne 


-40.42 ± 0.50 


12.7± 3.5 


[60] 


IN a 


— ijU.D it i.U 


Ll.i± i.D 


[61J 


24Mg 


-60.2 ± 1.2 


24.3 ± 1.6 


m 




-95.1 ± 2.0 


41.0 ± 4.0 


[62] 




-2402 ± 24 


441 ± 87 


[3]J 




-2483 ± 18 


326 ± 80 


[3]J 


i23Sn 


-2523 ± 18 


341 ± 72 


m 


205pb 


-5354 ± 61 


764 ± 165 


[63] 



The above potential is inserted into the KG equation ([21) to obtain a complex eigenvalue. 
The Coulomb potential due to the finite size charge distribution as well as the Uehling 
a{Za) vacuum polarization potential [s^ are also included. The strong interaction effects 
are the differences between these eigenvalues with and without the above potential (Eq. f[T^ ). 
respectively. 



B. Pionic atom data 



Experimental results on pionic atoms covering the whole of the periodic table have been 
pubhshed for a few decades, with improved accuracy over the years and with increased use of 
separated isotopes for targets. The most extensive analysis has been that of Konijn et al. 59l 
who analyzed 140 data points covering states from Is in to 4/ in ^^^Np. We note that a 



different definition of the strong interaction shift is used in Ref. [59j, namely, the difference 
between the complex binding energy for the full potential, including finite size Coulomb 
and vacuum polarization potential, and the binding energy for the point Coulomb potential. 
This is different from the conventional definition used here and, in any case, it breaks down 
for Is states when the charge of the nucleus is Z > 137/2. Such data are available now for 
the Is state in ^°^Pb. 

In the present work we are mostly concerned with global properties of the pion-nucleus 
interaction and its dependence on the nuclear density. This is effected by performing global 
fits to pionic atom data, handling together all relevant data. In order to avoid some distortion 
of the emerging picture we excluded from the analysis several deformed nuclei and also 
several very light nuclei, where the concept of an optical potential could be questionable. 
Moreover, we study extensively the dependence of strong-interaction effects on the neutron 
densities pn, using two-parameter Fermi distributions. For that reason we excluded from 
the data base also nuclei with Z < 8 where other densities such as a modified harmonic 
oscillator are more appropriate. 

The data used in the present work are summarized in the following four tables. The 
number of data points is 100, compared to 54 points used in our previous Review fi\. Ten 
of the points are due to the recently observed deeply bound states, as can be seen from the 
last four entries in Table [T] and the last line of Table [Tll These states are discussed below. 
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TABLE II: Data for 2p states in pionic atoms 





shift (keV) 


width (keV) 


Ref. 


24Mg 


0.129 




0.004 


0.0725 


±0.0018 


[62j 


26Mg 


0.126 




0.004 


0.0811 


±0.0019 


m 




0.286 




0.010 


0.192 


±0.009 


m 


30Si 


0.281 




0.010 


0.196 


±0.008 


m 


S 


0.615 




0.022 


0.430 


±0.021 


[65j 


40Ca 


1.941 




0.080 


1.590 


±0.023 


[6^ 


42Ca 


1.650 




0.080 


1.65 


±0.15 


mi 


^^Ca 


1.583 




0.080 


1.60 


±0.07 


[67] 


48Ca 


1.295 




0.115 


1.64 


±0.11 


[67] 


46xi 


2.490 




0.140 


2.39 


±0.15 


[611 


48xi 


2.290 




0.130 


2.62 


±0.15 


[63 


50Ti 


1.937 




0.110 


2.15 


±0.27 


[67] 


50cr 


3.560 




0.210 


4.3 


±0.4 


[6^ 




3.120 




0.190 


3.85 


±0.21 


[6^ 




2.827 




0.190 


3.84 


±0.29 


[68j 


Fe 


4.468 


± 


0.340 


6.87 


±0.21 


[6^ 


Ge 


5.5 




0.9 


18.5 ±2.5 


[69] 


As 


4.6 




0.9 


14.5 ±4.0 


[6^ 


Nb 


2.1 




3.0 


64 ± 8 


m 


Ru 


-27.6 ± 7.0 


77 ±24 


m 


205pb 


-835 ± 45 


321 ±61 


[63] 



C. Deeply bound pionic atom states 



As mentioned in Sec. [TTl the last decade has been dominated, for pionic atoms, by the 
observation of deeply bound states, which contributed to the revival of interest in pionic 
atoms in general and in the vriV in-medium interaction in particular. We therefore include 
here a brief outline of this topic. 

The term deeply bound pionic atoms refers to Is and 2p states in heavy pionic atoms 
which cannot be populated via the X-ray cascade process because upper levels such as the Sc? 
are already broadened by the strong interaction to the extent that the radiative yield becomes 
exceedingly small. The first to show that Is and 2p states in heavy pionic atoms, perhaps 
surprisingly, are so narrow as to make them well defined, were Friedman and Soff [75|. They 
calculated numerically binding energies and widths for pionic atom states well beyond the 
experimentally reachable region using optical potential parameters [t^ which reproduced 
very well experimental results over the entire periodic table. Figure H] is the original figure 
from 1985 (Ref. [tI]) showing the calculated binding energies {B) and widths (F) for Is 
states of pionic atoms as function of the atomic number Z. It is seen that up to the top end 
of the periodic table the widths of the states are relatively small. 

The explanation of that unexpected results was also given by Friedman and Soff 75|] in 
terms of the overlap of the atomic wavefunction with the nucleus. The atomic wavefunctions 
of these deeply bound states are excluded from the nucleus due to the repulsive s wave part 
of the potential such that their overlap with the imaginary part of the potential becomes 



16 



TABLE III: Data for 3d states in pionic atoms 





shift (keV) 


width (keV) 


Ref. 


Nb 


0.761 ±0.020 


0.402 ±0.016 


[62j 


Ru 


1.428 ±0.080 


0.75 ±0.08 


[62] 


Ag 


1.99 ±0.05 


1.43 ±0.04 


m 


Cd 


2.23 ±0.09 


1.65 ±0.07 


m 


Ba 


5.44 ±0.27 


4.3 ±0.9 


m 




7.05 ±0.30 


5.6 ±1.0 


m 


Ce 


7.23 ±0.30 


6.5 ±0.9 


m 


148Nd 


7.85 ±0.40 


8.8 ±1.2 




isoNd 


7.77 ±0.40 


9.2 ±1.1 




166gj. 


17.4 ±1.0 


19.7 ±0.9 


m 




16.3 ±0.5 


19.4 ±1.0 


m 


Pt 


28.8 ±1.9 


37 ±5 


m 


208pb 


32.25 ±2.40 


47 ±4 


m 


Bi 


30.5 ±3.0 


52 ±4 


[731 



TABLE IV: Data for 4/ states in pionic atoms 



shift (keV) width (keV) Ref. 

i^^Er 0.351 ±0.020 0.22 ±0.03 [74] 

Re 0.76 ±0.04 0.41 ±0.05 [73] 

Pt 1.13 ±0.04 0.59 ±0.05 [73] 

208pb 1.72 ±0.04 0.98 ±0.05 \J3\ 

Bi 1.83 ±0.06 1.24 ±0.09 [73] 

U 5.08 ±0.20 3.65 ±0.65 [65] 



very smalL In this context it is instructive to note that for a Schrodinger equation the width 
of a state is given exactly (i.e. not perturbatively) by the following expression 



~ f I / |2 . • v^"^J 

A slightly different expression is obtained for the KG equation. It is, therefore, easy to see 
that the reduced overlap of the atomic wavefunction with the nucleus, due to the repulsive 
real part of the s-wave term of the potential Eq. (fT5|) . provides a natural explanation of 
the numerically observed saturation of widths. Only for larger charge numbers, beyond the 
range of stable nuclei, the Coulomb attraction overcomes the s-wave repulsion, resulting in 
very large increase in the calculated widths. 

Three years later Toki and Yamazaki also concluded that deeply bound pionic atom levels 
would be sufficiently narrow and therefore could be observed, and in addition discussed 
experiments that could populate such states [77,]. After several unsuccessful attempts in 
various laboratories, it was realized that the key to success was the creation of a pion at 
rest (it appears that the first to suggest this approach, although in a different context, were 



Ericson and Kilian 78|). That requires 'recoil- free' kinematics, which for the (d,'^He) reaction 



on Pb creating a bound pion means a beam energy around 600 MeV. The experiment at 
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FIG. 4: Calculated binding energies and widths of Is states in pionic atoms (from Ref. j75l|). 



GSI used the fragment separator in order to achieve the required resolution and reduction of 
background and the first observation of the Is state in pionic atoms of ^°^Pb by Yamazaki 
et al. 3^ was a clear demonstration of the ability to study such states. Further experiments 
by the same group achieved improved accuracies for the Is and 2p levels in pionic atoms 



of ^^Pb using a ^"°Pb target [63] and studies of the Is state in isotopes of Sn [31[ followed 
soon after. Full details of this fascinating project can be found in a Review by Kienle and 
Yamazaki j^. 

To conclude this section, Fig. O illustrates the phenomenon of saturation of widths in 
pionic atoms by showing calculated radiation and total widths of 'circular' states (i.e. states 
with radial number n = Z + 1) in pionic atoms of Pb, using for the strong interaction 
the present-day best fit potential (see below). It is seen that for large I values the widths 
are essentially the radiative widths and X-ray transitions will be observed. As the /-value 
decreases the strong interaction width becomes significant and eventually, for low / values, 
the radiative transitions will be suppressed. However, the initial rise of the strong interaction 
width with decreasing / is saturated, resulting in the relatively narrow widths of deeply bound 
states. This is seen in Fig. showing the energy spectrum for circular states of pionic atoms 
of Pb. Comparing with Fig. it is clear that the energy levels are well-defined only thanks 
to the saturation of the widths. 
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FIG. 5: Calculated widths of 'circular' states in pionic atoms of Pb, showing radiative (dashed) 
and total widths (solid line). 
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FIG. 6: Calculated energy levels of 'circular' states in pionic atoms of Pb. The bars show the 
widths of the states. 
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D. Fits to pionic atom data 



1. General 



With nine parameters in the pion-nucleus optical potential it is not a straightforward 
task to get meaningful information on the potential from fits to the data. Moreover, 
there is the question of the proton pp and neutron p„ density distributions which are an 
essential ingredient of the potential. The proton density distributions are known quite well 
from electron scattering and muonic X-ray experiments [13], and can be obtained from 
the nuclear charge distributions by numerical unfolding of the finite size of the proton. In 
contrast, the neutron densities are not known to sufficient accuracy and their uncertainties 
must be considered when extracting parameters from fits to the data. 

Parameters of the pion-nucleus potential Eq. f|T^ had been obtained by performing fits 
to the data already in the late 1970s and early 1980s and with the additional accumulation 
of data they were found to possess good predictive power of yet unmeasured quantities, with 
little dependence on the details of the potential form used to fit the data. Two approaches 
have been made to the problem of large number of free parameters: (i) assumptions can be 
made in order to reduce the number of free parameters, such as introducing 'effective' density 
in the quadratic terms [79| thus making them linear in the density and avoiding correlations 
with the corresponding genuinely linear terms; (ii) performing global fits to very large sets 
of data where the effects of correlations between parameters are reduced. In most cases 
the dependence on neutron densities have not been studied and 'reasonable' densities have 
been assumed in the analysis. Only Garcia-Recio, Nieves and Oset 12| varied also neutron 
densities and obtained parameters for neutron distributions from global fits to pionic atoms 
data. Using three different versions for the pion-nucleus potential, they showed that the 
neutron rms radii were determined to good accuracy and were in reasonable agreement with 
values deduced from Hartree-Fock calculations. Below we present detailed studies of the 
dependence of derived potential parameters on neutron densities [8Q| , with emphasis on the 
's-wave anomaly' problem (see below). 

Addressing first the p-wave part of the potential (Eq. f|T6l) ). we note that due to its 
gradient nature it is expected to be effective only in the surface region of nuclei and that 
large medium or density-dependent effects are not expected beyond the 'trivial' LLEE cor- 
rection Eq. f[T7|) . As discussed in Ref. [l[, the dependence of on the parameters of the 
p-wave term is rather weak, with the free pion-nucleon scattering volumes of co=0.21m~^ and 
Ci=0.165m~^ being consistent with the best-fit values, when the 'classical' value of ^=1.0 is 
being used. We therefore adopt these three values and focus most of the attention on the s 
wave part of the potential. 

The role of the deeply bound states in the global picture of pion-nucleus interaction was 
an obvious question when the first such states were observed, and it was found [sii] that shifts 
and widths of deeply bound states agreed with predictions made with potentials based on 
'normal' states. That is fully understood by considering the overlaps between atomic and 
nuclear densities, which are very similar for normal and for deeply bound states, as can 
be seen e.g. from Fig. 17 of Ref. [3|. However, the experimental study of deeply bound 
states provides information on Is states throughout the periodic table and has enriched our 
undertanding of the pion-nucleus interaction at low energies. Studies of isotopes have the 
potential of providing additional valuable information. In what follows we combine the data 
on deeply bound states with data on normal states, as is seen in Tables [T] and [Tll 
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FIG. 7: Results of global fits to pionic atom data for different values of the neutron radius parameter 
7 of Eq. ([8]). Lower part - values of for 100 data points for three shapes of the neutron density 
Pn- Upper part - the corresponding isovector parameter h\ in comparison to its free -kN value 
(marked 'free'). 



2. The role of neutron densities 

As discussed in Sec. II C\ Vp, the rms radius for the proton density distribution, is con- 
sidered to be known and we therefore focus attention on values of the difference — Vp, 
using the linear dependence of r„ — Vp on (A^ — Z)/A given by Eq. ([H]). In order to allow for 
possible differences in the shape of the neutron distribution, the 'skin' and the 'halo' forms 
of Ref. [ll] were used, as well as their average, see Sec. II C[ 

Parameters of the potential Eq. (IT^ were determined by minimising the in fits to the 
combined data given in tables HI to II VI The linear p-wave parameters were held fixed at their 
respective free vrA^ values of co=0.21m~'^ and ci=0.165m~^ and the LLEE parameter was held 
fixed at ^=1. The quadratic p-wave parameters were found to be close to ReCo=0.01m~^ 
and ImCo=0.06m~^. The neutron radius parameter 6 was fixed at —0.035 fm and we scanned 
over the other radius parameter 7. 

Figure El shows results of global fits to the 100 data points for pionic atoms, as discussed 
above. The lower part shows values of function of the parameter 7, as defined by 

Eq. ([!])• It is seen that the quality of fit depends on the shape of the neutron distribution, 
where the 'skin' shape is definitely preferred. The existence of quite well-defined minima 
is gratifying. The upper part shows the corresponding values of the isovector parameter bi 
in comparison with its free ttN value (marked 'free'). It is clear that the resulting values 
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FIG. 8: Results of global fits for different values of the neutron radius parameter 7 of Eq. 
Lower part - the complex parameter Bq. Upper part - the isoscalar parameter 60 • 



of bi are almost independent of the shape of the neutron density used in the fit, depending 
mostly on the values of the rms radius r„, as represented by the parameter 7. From the 
figure one might conclude that for the best fit value of 7=1.6 fm the value of 61 turns out 
to be in reasonable agreement with the free nN value. 

Figure [8] shows in the upper part values of the isoscalar parameter bo- As noted above, 
this parameter is exceptionally small and the errors on the best-fit results are relatively 
large. The lower part of the figure shows the resulting values of the quadratic parameter 
Bq. The values of Im_Bo are determined to a very good accuracy and are independent of the 
neutron radius parameter 7 and of the shape of the distribution. The values of ReSo are not 
as accurate but they are clearly not zero, representing non-negligible repulsion in addition 
to the repulsion provided by the linear term. 

The rms radii of neutron distributions implied by the global fits to pionic atom data 
need some attention, particularly in view of the fairly strong dependence of the derived 
in medium values of the isovector s-wave parameter 61 on the rms radius assumed for p„. As 
mentioned above, values of the difference r„ — between rms radii obtained for medium- 
weight and heavy nuclei with 7=1.5 fm are too large by 0.05 to 0.1 fm in comparison with 



recent RHB calculations [21|, |22| . Moreover, a survey of different sources of information on 
^n — ^p suggests [16| that 7=0.9-1.0 fm is the proper value. Therefore the pionic atoms result 
of 7 > 1.6 fm is at odds with all we know about neutron density distributions in nuclei. 
An obvious next step is to try some finite-range (FR) modification of the otherwise zero 
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FIG. 9: Same as Fig. [7] but with finite range folding applied to the p-wave potential, see text. 



range pion-nucleus potential. This can be achieved by folding the densities with e.g. a 
Gaussian such that each density p in the potential is replaced by a folded one , see 
Eq. (fTUjl . It is found that when folded densities are used the becomes a monotonic 
increasing function of the range parameter f3. However, when folding is introduced separately 
into the s-wave and into the p-wave parts of the potential, then a minimum is obtained for 

when the rms radius of the finite range function is 0.9±0.1 fm, provided FR is applied 
only to the p-wave term. The same results are obtained when a Yukawa FR function is used 
instead of the Gaussian. Figure [9] shows results of global fits to pionic atom data when a FR 
folding with rms radius of 0.9 fm is applied to the p-wave part of the potential. Comparing 
with Fig. [7| we note that again the 'skin' shape yields the best fit but now it is obtained 
for 7=1.1 fm, which is an acceptable value. However, the difference between bi and its free 
ttN value is very large, representing a significant increased repulsion in the nuclear medium. 
Figure dU] is quite similar to Fig. [S] with non-zero values for Rei?o- Although there is no 
acceptable theory for the empirical parameter Bq, it is believed that the absolute value of 
its real part should be smaller than the imaginary part (82], 8^], which is not the case here. 

The present results mean that the in-medium s-wave repulsion, as obtained from global 
fits to pionic atom data, is significantly enhanced compared to expectations, partly via the 
extra repulsion provided by bi and partly by the dispersive ReBo being more repulsive than 
expected. The sum of these two effects is the well-known pionic atoms 'anomaly'. This is 
the most bias- free way of presenting the 'anomaly', or 'anomalous s-wave repulsion'. 
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FIG. 10: Same as Fig. [8] but with finite range folding applied to the p-wave potential, see text. 
3. The s-wave anomaly and the issue of chiral symmetry restoration 

The so-called s-wave anomaly or the extra repulsion observed in fits to pionic atom data 
had been known for a very long time pLj]. The previous section presented a state-of-the-art 
summary of this topic, based on global fits to 100 data points, respecting our knowledge 
on the r„ — rms radii difference and otherwise without additional assumptions. Global 
approaches always yielded the result that the extra repulsion observed in the phenomenolog- 
ical potential is due to two sources, namely, an enhanced bi coefficient and an unexpectedly 
large repulsion of the dispersive quadratic term, albeit with large uncertainty. An alterna- 
tive approach, mainly due to Yamazaki and co-workers, has been to handle restricted data 
sets and to reduce the number of parameters in the potential by making some assumptions. 
Basically they used the approach of Seki and Masutani j79| where due to the correlations 
between and Rc-Bq and the assumption of an average or an effective density, the two terms 
plus the isoscalar double-scattering contribution of Eq. (^1]) are lumped together, resulting 
in a single effective isoscalar real part linear in the density. Using the deeply bound Is 
states in Sn isotopes (see Table [T]) together with Is states in ^^O, ^°Ne and ^^Si, a total of 
12 points, they obtain bi = —0.1149 ± 0.0074m~^ (see Ref. [1] for full details.) This may be 
compared with the value of bi that corresponds to the minimum of in Fig. [9] for the 'skin' 
shape of pn, namely bi = —0.109 ± O.OOSm^^. The agreement between the two results is 
very good, and a significant discrepancy is thus established beyond any doubt with respect 
to the free nN value b{ = —0.0864 ± O.OOlOm"^ derived from the preliminary results of the 
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PSI measurement of the vr H Is level shift and width 8J]. Note, however, the conceptual 
differences between the two methods: 

• The Yamazaki et al. approach makes particular assumptions on the neutron densities 
of only few nuclides for Is atomic states. In our global approach, where it is found that 
the major effect of the unknown neutron densities is through r„ — r^, the (A^ — Z)/A 
dependence of r„ — is an average over 36 nuclides spanning a full range of atomic 
states (from Is to 4/). 

• Due to the use of 'effective density' as mentioned above, and N = Z nuclei dominated 
by the isoscalar vrA^ interaction, it is somewhat ambiguous in the Yamazaki et al. 
approach to determine 60 and Re-Bo independently of each other. Thus, assuming 
60 = 0, they obtain Re^o = -0.033 ± 0.012m-^ [Eq. (77) of Ref. Q]. Our data base 
includes 100 rather than 12 points, where all are significant in determining parameters 
of the s-wave part of the potential since this part contributes large fractions of the 
strong interaction effects also for states with / > 0. Moreover, no assumptions are 
made regarding the terms nonlinear in density. In addition figure (TU] shows that our 
deduced value of bo for the 'skin' shape for p„ is essentially in agreement with the free 
vrA^ value feg = +0.0068 ± 0.0031m~^, deduced from the preliminary PSI results for 
the TT^H Is level shift and width |84|]- If this value were assumed in Ref. [3|, then they 
would have derived a value of Rei?o = —0.064 ± 0.012m~^ agreeing within error bars 
with our values as depicted in Fig. [TOl 

A systematic study of the uncertainties in parameters of the potential and their depen- 
dence on the size of the data base and on assumptions made and constraints imposed in 



the analysis can be found in Ref. [17|. In what follows we focus attention on the in- medium 
values of 61. It is shown below that eventually the problem with Rc-Bq is also solved. 

The renewed interest in recent years in the 'anomalous' s-wave repulsion in the pion- 
nucleus interaction at threshold, as found in phenomenological analyses of strong interac- 
tion effects in pionic atoms, is partly due to the suggestion by Weise jsl, 8^ that such 



enhancement could be expected, at least in the isovector channel, to result from a chirally 
motivated approach where the pion decay constant becomes effectively density dependent in 
the nuclear medium. Since bi in free-space is well approximated in lowest chiral-expansion 
order by the Tomozawa- Weinberg (TW) expression [Sj 



bi = = -0.08m;^ , (24) 

with /x^jv the pion-nucleon reduced mass, then it may be argued that bi will be modified in 
pionic atoms if the free-space pion decay constant = 92.4 MeV is modified in the medium. 
QCD coupled with PCAC relates to the quark condensate < qq >: 

/^M = liiZ£^l_4^L, (25) 

<qq>o mlf^ 

where a is the nN sigma term and where the last step provides the leading term in a density 
expansion of the quark condensate [89| assuming that the charge-averaged pion mass does 
not change in the medium.^ Thus /^r is expected to decrease in the nuclear medium by 



2 



This holds for the temporal version of /tt, corresponding to the vacuum-to-pion matrix element of the 



time component of the axial current, c.f. Refs. |9d.l91|. 
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FIG. 11: Results of global fits for different values of the neutron radius parameter 7 (Eq. ([8])) 
with bi{p) given by Eq. (j26p . Lower part - values of for 100 data points for three shapes of the 
neutron density pn. Upper part - the corresponding isovector parameter bi in comparison to its 
free vrA'^ value (marked 'free'). 



about 20% at nuclear-matter density po — 0.16 fm ^, using a ^ 50 MeV 92]. The form of 



Eq. (HM suggests then a density-dependent isovector amplitude such that bi becomes 

1 - ap/mif^ 1 - 2.3p 

for (7=50 MeV and with p in units of fm~^, resulting in an increase of bi in the nuclear 
medium by about 60% at nuclear-matter density po — 0.16 fm~^. This ansatz 85, 8^ was 



applied at almost the same time in two different analyses of pionic atom data. Kienle and 
Yamazaki ^] outlined a method for extracting bi from analyses of very limited data sets, 
which is essentially the method described above, but using a fixed average value of the 
density-dependent &i(p) of Eq. (1261) . Friedman [9^ presented results of global analyses, as 
outlined above but with 60 data points compared to the present 100 points, using explicitly 
Eq. fl26p for the density dependent &i(p). It was shown that indeed most of the difference 
between the derived bi and its free ttN value disappeared when the above density dependence 
was included. 

Figure [TT] shows results of global fits to the 100 data points for pionic atoms, with 61 (p) 
given by Eq. fl26l) . The parameter bi stands here for 61 (p = 0). It is seen that for the 
lowest minimum of x^? i-e. the minimum of the curve obtained for the 'skin' shape for the 
neutron density, bi is much closer to the free 71 N value than it was in Fig. [9] where a fixed bi 
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FIG. 12: Results of global fits for different values of the neutron radius parameter 7 (Eq. ([8])) 
with bi{p) given by Eq. ([26]) . Lower part - values of the complex parameter Bq. Upper part - the 
resulting isoscalar parameter bo. 





was assumed, and they are almost in agreement. Figure [T^ shows that bo now is consistent 
with zero and that Re-Bo is much less repulsive than before and is acceptable being almost 
consistent with zero. 

The renormalization of in dense matter, as given by Eq. fl25l) and leading to the 
related renormalization of the isovector amplitude bi in Eq. fl26|) . may also be derived under 
simplifying assumptions by solving the KG equation in infinite nuclear matter of protons 
and neutrons for a pion chiral s-wave polarization operator Ils{E) near threshold [oS, looj. 
In chiral-perturbation theory, in the limit of vanishing mj^/M and for zero momentum pions, 
q = 0, the leading terms to order E/lAirfT^) of the vrA^ s-wave amplitudes give rise to the 
following form of the s-wave pion polarization operator: 

n,(E) = 2EV}$iE) ^ r,-f^(p„ - p,) - ^^^|^(Pp + Pn) , (27) 

where E is the pion energy including its rest mass and = +1,0,-1 for tt", vr'', tt"*", 
respectively, in the notation of Eq. ([2]). The first, TW isovector term is the dominant one 
near threshold, providing repulsion for vr" for all nuclei with neutron excess. The second, 
isoscalar term is nearly zero at threshold, so it is reasonable to fit (3 by requiring f3 = (j/m^. 
The polarization operator 'n.s{E) satisfies a KG equation (for q = 0): 

E^-m^ -n,(E) = (28) 
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which by inserting Eq. (127|) becomes 



)(E'-ml)-Tz^{pn-p,) = 0, (29) 



TV J TV J TV 

with p = Pp + Pn- When recast into the form 

n,(E) = E^-ml = T.^ip,, - p,) , (30) 
with the effective density-dependent pion decay constant /* defined by 

/^^/.^(l-^), (31) 

it is seen to be equivalent, at threshold, to Eqs. (12411261) with the renormalized effective 
pion decay constant /*. Equation (!30!) strictly speaking is satisfied for one specific value of 
energy ('self energy') near threshold, E ^ m^r. Approximating then Ils{E) by the right-hand 
side at = m^, the appearance of /* instead of /^r may be attributed to a wave- function 
renormalization effect josl . 96| . 



Switching from infinite nuclear matter considerations to actual pionic-atom calculations, 
the finite-size Coulomb potential Vc needs to be introduced properly into a KG equation in 
which Ils{E) serves as a given input, not as an entity to solve Eq. fl28l) for. Schematically, 
Eq. fl28|) is replaced by 

(E - V,f -ml- Iis{E - \/,)j ^ = , (32) 



where the Coulomb potential Vc enters via the minimal substitution requirement [9j 
H-siE) I{s{E — Vi). When the chiral version of the s-wave pion polarization operator 
Eq. was used in a KG equation of the type of Eq. ([32D, it was found in the global fits 
to pionic atom data reported in Ref. [ol] that a large over- correct ion of hi occurred, to a 
value hi = —0.068 ± 0.004m~^, significantly less repulsive than 6{. In this case, it is the 
combined effect of the isoscalar and the isovector amplitudes that is responsible for over- 
shooting Therefore, although the effect of including the energy dependence of the q = 
chiral amplitudes goes in the desired direction, it does not provide a quantitative resolution 
of the s-wave anomaly problem. Different results, and conclusions, are obtained when fits 
are made to partial data sets which may not carry sufficient statistical significance to decide 



one way or another on this issue [95|, |96 



Friedman and Gal [98| discussed also an alternative procedure in which the empirical on- 
shell energy dependent ttN amplitudes are used for implementing the minimal substitution 
requirement E ^ E — Vc- As pointed out by Ericson [99| the on-shell approximation 
follows naturally for strongly repulsive short-range NN correlations from the Agassi-Gal 



theorem {100| for scattering off non-overlapping nucleons. The energy dependence of the 
'empirical' amplitudes is weaker than that of the 'chiral' ones. For hi{E) there is hardly any 
energy dependence, in contrast to the energy dependence of the 'chiral' hi{E) which by itself 
would suffice to recover 6(, the free-space value of hi, in the pionic-atom global fits. For 
bo{E) the empirical energy dependence is only about 60% of the 'chiral' effect. When this 
weaker empirical energy dependence of ho was used, the resulting hi was still too repulsive 
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FIG. 13: Results of global fits for different values of the neutron radius parameter 7 (Eq. ([8])) with 
bi{p) given by Eq. ([26]) and the empirical energy dependence of the free bo{E) included. Lower 
part - values of for 100 data points for the 'skin' shape of the neutron density pn- Upper part - 
the corresponding isovector parameter bi in comparison to its free vrA^ value (marked 'free'). 

in comparison with but closer to b{ than the 61 resulting in the conventional fixed bi 
model. 

Figure [13] is similar to Fig. [11] but with the energy dependence of the free vrA^ 60 included 
in addition to applying the density-dependent renormalization of Eq. f[26l) for hi. Comparing 
with Fig. [11] the values of bi have shifted now and at the minimum of the agreement with 
the free vrA^ value is perfect. The resulting 60 (not shown) of 60 = — 0.009 ±0.007m~^ is close 
to the free value and KeBo (not shown) is essentially zero at — 0.005 ±0.035m~'^. The chiral- 
motivated isovector 61 (p) Eq. fl26|) together with the empirical energy dependence therefore 
provides all the required extra repulsion, with the isoscalar bo being essentially zero and 
with no significant dispersive Rc-Bq term needed. 



4- Radial sensitivity of pionic atoms 

The chiral-motivated dependence of the isovector term bi, as given by Eq. fl26p . almost 
completely removed the 'anomaly' observed when a fixed value was used for bi. Not only 
with Eq. fl26p does now one get for 61 (0) the free vrA^ value, but also the fit value of 60 is 
very close to its free value and Rei?o is zero. This is obtained without the need to introduce 
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'effective' density or to make any assumptions regarding ReBo- It is concluded, therefore, 
that Eq. (126|) is a fair representation of the medium-modification of the isovector term of 
g(r), Eq. (ITSl) . It is instructive to examine further the radial sensitivity of this term. 

The radial sensitivity of exotic atom data was addressed before [l[ with the help of a 
'notch test', introducing a local perturbation into the potential and studying the changes in 
the fit to the data as function of position of the perturbation. The results gave at least a 
semi-quantitative information on what are the radial regions which are being probed by the 
various types of exotic atoms. However, the radial extent of the pert urbation was somewhat 



arbitrary and only very recently that approach was extended |lOll | into a mathematically 
well-defined limit. 

In order to study the radial sensitivity of global fits to exotic atom data, it is necessary 
to define the radial position parameter globally using as reference e.g. the known charge 
distribution for each nuclear species in the data base. The radial position r is then defined 
as r = i?c + V(^c, where Rc and Oc are the radius and diffuseness parameters, respectively, of 
a 2pF charge distribution T^. In that way rj becomes the relevant radial parameter when 
handling together data for several nuclear species along the periodic table. The value of 
is regarded now as a functional of a global optical po tential V( ri), i.e. = X^[V{v)]y where 
the parameter 77 is a continuous variable. It leads to (lOll . 102l | 



where 

5V{ri') a~>oev^o ev 

is the functional derivative (FD) of The notation 5a{f]—ri') stands for an approximated 

(5-function and ey is a change in the potential. From Eq. ( l33l) it is seen that the FD 
determines the effect of a local change in the optical potential on x^. Conversely it can 
be said that the optical potential sensitivity to the experimental data is determined by the 
magnitude of the FD. Calculation of the FD may be carried out by multiplying the best fit 
potential by a factor 

/ = l + e5.(r/-V) (35) 
using a normalized Gaussian with a range parameter a for the smeared 5-function, 



5.(r/-V) = ^e-(''-''')'/2'^^ (36) 



(37) 



For finite values of e and a the FD can then be approximated by 

5x^\vm 1 x'[vmi + e5^{r^ - j))] - x^[vm 

5V{7]') ^ V{7]') e 



The parameter e is used for a fractional change in the potential and the limit e — is 
obtained numerically for several values of a and then extrapolated to cr = 0. 

Figure [m shows the FDs with respect to relative changes in the real part and with respect 
to relative changes in the full complex s-wave part of the best-fit pion-nucleus potential 
where the chiral 61 (p) is assumed. Calculations of the FD with respect to the imaginary 
part of the potential show additivity of the FDs, hence the difference between the FD for 
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FIG. 14: Functional derivatives of the best fit ^ with respect to relative changes in the real part 
(solid curve) and with respect to relative changes in the full complex (dashed) s-wave pion nucleus 
potential Eq. with 6i(p) of Eq. ([26]) 



the full complex potential and for the real part is the FD with respect to relative changes 
in the imaginary part of the s-wave part of the potential. It is immediately clear from 
the figure that both parts of the potential have quite similar importance in determining 
strong interaction effects in pionic atoms. Turning to radial sensitivity, recall that 77 = —2.2 
corresponds to 90% of the central density of the nuclear charge and ri = 2.2 corresponds to 
10% of that density. It is therefore clear that pionic atom data are sensitive to the s-wave 
part of the potential over a region where the densities change between free space to the 
full central density. It is no wonder that using a fixed value for bi in fits to data led to 
'abnormal' values for this parameter. With the success of the chiral &i(p) there is no reason 
for employing approximations such as linearization, by using an effective density. 



E. Pion elastic scattering 



The applicability of the Kisslinger potential [52(] and its Ericson-Ericson extension [55 



(Eqs. (11511221) ) to elastic scatte ring of pions by nuclei had been demonstrated early in the 



days of the 'pion factories' [79|, |103| . mostly at energies higher than 40-50 MeV. Following 
the success with the 'chiral' density dependence of bi in the subthreshold regime of pionic 
atoms, it is interesting to study the behavior of the pion-nucleus potential across threshold 
into the scattering regime and to examine whether the pionic atom 'anomaly' is observed also 
above threshold when using a fixed value for the parameter bi. Of particular interest is the 
question of whether the density dependence of the isovector amplitude as given by Eq. fl26|) 
is required by the scattering data. In the scattering scenario, unlike in the atomic case, one 
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FIG. 15: Experimental results for elastic scattering of vr^ compared with calculations. 



can study both charge states of the pion, thus increasing sensitivities to isovector effects and 
to the energy dependent effects due to the Coulomb interaction. Looking for earlier suitable 
data for elastic scattering, it was somewhat surprising to realize that at kinetic energies well 
below 50 MeV there seemed to be only one set of high quality data available for both c harge 



states of the pion obtained in the same experiment, namely, the data of Wright et al. [104 



for 19.5 MeV pions on calcium, with predominantly the N = Z isotope Ca. For that 
reason precision measurements of elastic scattering of 21.5 MeV 7r+ and tt" by several nuclei 
were performed very recently at PSI (32I . 33 1 with the aim of analyzing the results in terms 
of the same effects as found in pionic atoms. The experiment was dedicated to studying 
the elastic scattering of both pion charge states and special emphasis was placed on the 
absolute normalization of the cross sections, which was based on the parallel measurements 
of Coulomb scattering of muons. 

Figures [15] and [16] show comparisons between experimental differential cross sections for 
the elastic scattering of 21.5 MeV vr"*" and vr" by Si, Ca, Ni and Zr and calculations made 
with the best-fit EW potential, defined below. The fit to the data was made to all eight 
angular distributions put together, see Ref. [3^ for details. Table |V] summarizes values of 
bi obtained for the various models, with C indicating a conventional (fixed bi) potential, W 
indicating the use of the chiral-motivated of Eq. f[26l) and EW stands for using &i(p) 
and including also the empirical energy dependence of bo- It is evident that with the C 
potential using a fixed value for bi the fit to elastic scattering results yields too repulsive a 
value for bi in comparison with the corresponding value for the free ttN interaction, much 
the same as is the case with pionic atoms. When the chiral motivated 61 (p) of Eq. (12^ 
is used is reduced significantly and the resulting bi is essentially in agreement with the 
free value. Adding the empirical energy dependence of 60 brings the two values into full 
agreement, much the same as is the case with pionic atoms. As with pionic atoms, the 
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FIG. 16: Experimental results for elastic scattering of vr^ compared with calculations. 



TABLE V: Values of bi from fits to elastic scattering of 21.5 MeV vr^ by Si, Ca, Ni and Zr. 
The free ttN value is b( = -0.0864 ± O.OOlOm-^ 



model 


C 


W 


EW 


6i(m-i) 


-0.114 ±0.006 


-0.081 ± 0.005 


-0.083 ± 0.005 


for 72 points 


134 


88 


88 



chiral energy dependence in Eq. (!27j) caused large over-correction of the resulting bi which 
turns out to be significantly less repulsive than its free nN value. 

It is therefore concluded that the extra-repulsion or s-wave anomaly is observed also in 
the scattering regime and that the medium-modification of the isovector term of the local 
part of the pion-nucleus potential is of the same nature both below and just above threshold. 

F. Conclusions 

Pionic atoms form the oldest type of exotic atom of a strongly-interacting particle and 
both experiment and theory had achieved maturity already in the 1990s, with the nagging 
problem of anomalously enhanced repulsion in the s-wave part of the potential essentially 
unresolved. The discrepancy between the in-medium bi and its free irN value became clearer 
with the ever improving accuracy of the experimental results on pionic hydrogen. In contrast, 
it seemed that experimental X-ray spectroscopy of medium-weight and heavy pionic atoms 
reached a dead end without new experiments in the last 20 years. The predicted existence 
of well-defined Is and 2p states in heavy pionic atoms, due to saturation of widths caused by 
the s-wave repulsion, prompted the pioneering (d,^He) experiments which supplied strong 
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interaction data unreached otherwise. Although the overall picture of pionic atom potentials 
has not changed due to these new data, they gave added impetus to the study of medium 
effects in the pion-nucleus interaction at threshold. The issue of the enhanced repulsion in 
the s-wave part of the potential appears to have been resolved with the chiral-motivated 
61 (p) of Eq. fl26l) which was shown to apply also at 21 MeV thanks to a recent dedicated 
experiment on elastic scattering of vr^ by nuclei. It may be concluded that the dominant 
effect of the nuclear environment on the real part of the pion-nucleus potentials close to 
threshold is given by Eq. (|26|) . Note, however, that the absorption terms quadratic in 
density remain largely phenomenological. 



IV. K NUCLEAR PHYSICS 



Preview 



The ^-nucleus interaction near threshold is strongly attractive and absorptive as sug- 
gested by fits to the strong- interaction shifts and widths of i^'-atom levels [!]• Global fits 
yield 'deep' densi ty deperident optical p otentials with nuclear-matter depth ReV^(po) ~ 
-(150-200) MeV [lOll . [lOSl . \wK llOTl . llOq , whereas in the other extreme case several studies 
that fit the low-energy K~p reaction data, including the / = quasibound state A (1405) 
as input for constructing self consistently density dependent optica l potentia l s, ob tain rel- 



lOgl . ITTol . [m|. For a 



atively 'shallow' potentials with ReV^(po) ~ -(40-60) MeV [2_, ^ 

recent update of these early calculations, see Ref. II2I ]. An example of a chirally motivated 
coupled-channel fit to the low-energy K~p cross sections is shown in Fig. [17] from Ref. [24 



This calculation is based on the free-space chiral coupled-channel amplitudes used i n th e 



in-medium calculations of Waas et al. |26] following the earlier work of Kaiser et al. [113 
By imposing self consistency in its nuclear part (solid lines) the calculation jl^ leads to a 
weakly density dependent shallow ^-nucleus potential in terms of the effective scattering 
length aeff(p) of Fig. [3] in Sec. [11 As is shown below, 'shallow' potentials are substantially 
inferior to 'deep' ones in comprehensive fits to i^'-atom data. The issue of depth of KeVj^ 
is reviewed below and the implications of a 'deep' potential for the existence and properties 
of X-nucleus quasibound states are discussed. Since the two-body KN interaction provides 
a starting point in many theoretical works for constructing the K nuclear optical potential 
Vk, we start with a brief review of the KN data available near the K~p threshold. 



B. The K p interaction near threshold 



The K~p data at low energies provide a good experimental base upon which models for 
the strong interactions of the KN system have been developed. Near threshold the coupling 
to the open vrE and vrA channels is extremely important, as may be judged from the size 
of the K~p reaction cross sections, particularly K~p vr+S", with respect to the K~p 
elastic cross sections shown in Fig. [T71 Theoretical models often include also the closed 



7] A, rjT,, K'E channels [114| . Other threshold constraints are provided by the accurately 



determined threshold branching ratios [ll5l . III6 
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FIG. 17: Calculations from Ref. 2j] of cross sections for K~p scattering and reactions. The dashed 
lines show free-space chiral- model coupled-channel calculations using amplitudes from Ref. 26l |. 
The solid lines show chiral-model coupled-channel calculations using slightly varied parameters in 
order to fit also the K~-atom data for a (shallow) optical potential calculated self consistently. 
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Additional sources of experimental information are the vrS invariant mass spectrum in the 
/ = channel from various reactions, and the K~p scattering length deduced from the 
recent measurements at KEK [3^ and at Frascati [3^1 (the DEAR collaboration) using a 



Deser-based formula 117 



(41) 
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FIG. 18: Calculated cross sections for K~p vr^E^'^ multiplied by 4|q^ ^'|y^ and continued below 
the K~p threshold (vertical line), for three chiral coupled-channel fits to the K~p low-energy data. 
The fit shown b y the solid (dashed) lines excludes (includes) the DEAR value for a^-p- Figure 



taken from Ref. 119 ]. 



where a is the fine-structure constant. The value of ax-p derived from the DEAR measure- 
ment using this expression, ax-p = (—0.45 ± 0.09) -|- i(0.27 ± 0.12) fm, appears inconsistent 
with most comprehensi ve fits to the l ow-e nergy K~p scattering and reaction data, as dis- 
cussed by Bora soy e t al. 118 , liol . 120, 121 |. A dissenting view, however, is reviewed recently 
by Oiler et al. fl2^ . 

The K~p elastic, charge-exchange and reaction cross sections shown in Fig. [17] refer to 
energies above the K~p threshold at ^/s = 1432 MeV. However, by developing potential 
models, or limiting the discussion to phenomenological K-matrix analyses, KN amplitudes 
are obtained that allow for analytic continuation into the nonphysical region below threshold. 
Using a K-matrix analysis, this was the way Dalitz and Tuan predicted the existence of 



the A(1405) vrS, 1 = reso nance in 1959 |123l |. Recent examples from a coupled- channel 



potential model calculation [ll9|| based on low-energy chiral expansion are shown in the 
next three figures. Figure [T8] depicts the quantity 4:\q^^P\y/saK-p^TT^j:±{s) obtained by 
continuing the K~p —>■ ttS amplitudes below the K~p threshold. One sees clearly the 
resonant behavior of the extrapolated cross sections due to the A(1405). This resonance 
is also seen in Fig. [19] where the real and imaginary parts of the K~p elastic scattering 
amplitude, continued analytically below the K~p threshold, are plotted. The discrepancy 
with Ima^-p deduced from the DEAR measurement 35|, as given above, is clearly seen. In 
contrast to amplitudes which allow for the / = channel and thus exhibit a resonance 
effect due to the A(1405), the purely 1 = 1 K~n amplitude does not show such effects below 
threshold, as seen in Fig. [20] where the real and imaginary parts of the K~n elastic scattering 
amplitude, within the same coupled-channel model, are shown. The model dependence of 
the K~n elastic scattering amplitude, as given by the three different curves, is considerably 
weaker than the model dependence of amplitudes in which the A (1405) resonance enters, 
e.g. the K~p tt^E^ amplitudes related to the cross sections shown in Fig. [TS] 
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FIG. 19: Real and imaginary parts of the K~p forward elastic scattering amplitude, fitted within 
a NLO chiral SU(3) coupled-channel approach to K~p scattering and reaction data. The line 
denoted WT is the (real) LO Tomozawa- Weinberg K~p driving-term amplitude. The DEAR 



measurement 
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35| value for a^-p is shown with error bars. Figure taken from Ref. 124i |. based on 
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FIG. 20: Real (left) and imaginary part (right) of the K^n forward elastic scattering amplitude, 
continued below the K~n threshold (vertical line). The various lines correspond to different in- 
teraction inodels, the dashed line standing for the Tomozawa- Weinberg interaction. Figure taken 
from Ref. [ui]. 



C. nucleus potentials 



The gross features of low-energy KN physics, as demonstrated in the previous section by 
chiral coupled-channel fits to the low-energy K~p scattering and reaction data, are encap- 
sulated in the leading-order Tomozawa- Weinberg (TW) vector term of the chiral effective 
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Lagrangian 125| . The Bom approximation for the i^-nwc/ear optical potential due to 
the TW interaction term yields then a sizable attraction: 



V, 



K 



P 



-55 — MeV 
Po 



(42) 



for po = 0.16 fm^^, where /tt ~ 93 MeV is the pseudoscalar meson decay constant. Iter- 
ating the TW term plus next-to-leading-order terms, within an in-medium coupled-channel 
approach constrained by the KN — vrE — vrA data near the KN threshold, roughly doubles 
this ^-nucleus attraction. It is found (e.g. Ref. |2^) that the A(1405) quickly dissolves in 
the nuclear medium at low density, so that the repulsive free-space scattering length aK-p, 
as function of p, becomes attractive well below po- Since the purely 1=1 attractive scat- 
tering length QK-n is only weakly density dependent, the resulting in-medium KN isoscalar 
scattering length 6o(p) = \{.0'K-p{.p) + ciK~n{p)) translates into a strongly attractive V^'- 



VRir) 



27r 



Pkn 



Mp) Pir) , ReVKipo 



-110 MeV 



(43) 



This in-medium KN isoscalar scattering length &o(p) is the same one as the effective scatter- 
ing length OefT plotted as function of the density p in Fig. [3] in Sec. II D[ It is useful to record 
that &o(po) ~ 1 fm corresponds to Vr{po) ~ —100 MeV. However, when Vr is calculated 
self consistently, namely by including Vr in the propagator Go used in the Lippmann- 
Schwinger equation determining 



ReVj^(po 



-(40-60) MeV |24j. Il09l. IllOl. Illl 



bnip), as demonstrated by Eq. (fT3|) in Sec. II Dt one obtains 



The main reason for this weakening of Vr, 
approximately going back to that calculated using Eq. is the strong absorptive effect 



which Vr exerts within Gq to suppress the higher Born terms of the KN TW potential. 
Additional considerations for estimating Vr are listed below. 



QCD sum-rule estimates 126l | for vector (v) and scalar (s) self-energies: 
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- 30 MeV , 



(44) 
(45) 



where is the strange-quark (current) mass. The factor 1/2 in Eq. (jUj) is due to the 
one nonstrange antiquark q in the K meson out of two possible, and the minus sign is 
due to G-parity going from q to q. This rough estimate gives then V^(po) ~ —130 MeV. 



The QCD sum-rule approach essentially refines the mean- field argument 127 . 128l 



Vk{Po) 



3 (Ss(iV) 



170 MeV 



(46) 



where the factor 1/3 is again due to the one nonstrange antiquark in the K meson, 
but here with respect to the three nonstrange quarks of the nucleon. 

The ratio of K~ / K^ production cross sections in nucleus- nucleus and proton- nucleus 
collisi ons near threshold, measured by the Kaon Spectrometer (KaoS) collabora- 
tion 129l | at SIS, GSI, yields an estimate V^(po) ~ —80 MeV by relying on 
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TABLE VI: Calculated {K^^^p,-ir) rates on ^^C per stopped K (in units of 10 ^) for Ipjy — > ls\ 
capture into the 1~ ground-state configurations in an d P B, for various fitted optical potentials 



ordered according to their depth. Table taken from Ref. 24\ . 



final 


chiral 


effective 


fixed 


DD 


Rev X 103 


12p 
A ^ 


0.231 
0.119 


0.169 
0.087 


0.089 
0.046 


0.063 
0.032 


0.98 ±0.12 [130] 
0.28 ±0.08 [131] 



Boltzmann-Uehling-Uhlenbeck (BUU) transport calculations normalized to the value 
Vk{po) ~ ±25 MeV. Since KNN YN absorption processes apparently were dis- 
regarded in these calculations, a deeper may follow once nonmesonic absorption 
processes are included. 

Capture rates of (-ft'stop! ^) reactions to specific A hypernuclear states provide a sensi- 
tive measure for the strength of the K~ optical potential at threshold. A very strong 
potential, as discussed in Sec. II VE] generates nwc/ear quasibound states which due 
to orthogonality with the K~ atomic states force the wavefunctions of the latter to 
oscillate and become suppressed within the nuclear volume, and thus to reduce their 
effectiveness in the calculation of the (-ftTstop; ^r) transition matrix element. The results 
of DWIA calculations jl^ for capture on ^^C into the 1~ ground state configurations 
of ji^^C and j^^B are shown in Table IVll using several fitted K~ optical potentials which 
are ordered according to their depth. It is clearly seen that the deeper the K~ optical 
potential is, the lower the calculated rate becomes. Unfortunately, all the calc u lated 



rates shown in Table |Vl] are much lower than the measured values -Rexp |l3Cll . 131 
making it difficult to draw definitive conclusions. Furthermore, the experimental cap- 
ture rates depart strongly from the ratio 2 : 1 expected from charge independence for 
the ratio R{fC)/R{fB) for the rates i?exp- 



D. Fits to K -atom data 



The K~-aXom. data used in global fits [1] span a range of nuclei from ^Li to ^^^U, with 
65 level-shifts, widths and transition yields data points. It was shown already in the mid 
1990s [1] that although a reasonably good fit to the data is obtained with the generic tp 
potential of Eq. ([3]) with an effective complex parameter ho corresponding to attraction, 
greatly improved fits are obtained with a density-dependent potential, where the fixed ho 
is replaced by ho ± Bo[p{r-) j po]"^ , with ho, Bo and a > determined by fits to the data. 
Fitted potentials of this kind are marked DD. This parameterization offers the advantage 
of fixing ho at its (repulsive) free-space value in order to respect the low-density limit, while 
relegating the expected in-medium attraction to the Bo term which goes with a higher 
power of the density. As mentioned in Sec. IIV A[ the tp best-fit potentials have real parts 
which are less than 100 MeV deep for medium-weight and heavy nuclei. The corresponding 
density-dependent potentials are more attractive, 150-200 MeV deep, hence the 'shallow' vs 
'deep' terminology. Chirally inspired approaches that fit the low-energy K~p reaction data 
predict attractive potentials of depths ~(110-120) MeV [26| and, additionally, by requiring 
self consistency in the construction of the optical potential, lead to yet shallower potentials 
with ReV^(po) ^ 



(40-60) MeV |24j . Ill0| . Recent experimental reports on candidates for 



39 




FIG. 21: Left: real part of the iC- Ni potential obtained in a global fit to K -atom data using the 



model-independent FB technique 10l[ |. in comparison with other best-fit potentials and values 
in parentheses. Right: con tributions to the of atomic shifts for the deep density-dependent 



potential F from Ref. 108t | and for the shallow tp potential. 



i^-nuclear deeply bound states in the range of binding 5^ ~ 100-200 MeV 0, H M, 113 
again highlighted the question of how much attractive the i^-nucleus interaction is below 
the KN threshold. Therefore the motivation for re-analysis of a comprehensive set of kaonic 
atom data is two fold. First is the question of 'deep' vs. 'shallow' real ^-nucleus potential, 
in light of recent possible experimental evidence for the existence of deeply bound kaonic 
states whose binding energies exceed the depth of the shallow type of potential obtained from 
fits to kaonic atom data. However, if the deep variety of potential is confirmed, then the 
dependence of the KN interaction on the nuclear density becomes of prime concern. This 
densi ty d ependence is the second point which motivated the re-analysis of kaonic atoms 



data [1081 ] although there have not been any new experimental results on strong interaction 



effects in kaonic atoms since the early 1990s. 

The departure of the optical potential from the fixed-t tp approach was introduced 108 
with the help of a geometrical model, where one loosely defines in coordinate space an 
'internal' region and an 'external' region by using the multiplicative functions F{r) in the 
former and [1 — F{r)] in the latter, where F{r) is defined as 

F(r) = (47) 

with X = {r—Rx)/ax. Then clearly F(r) 1 for (i?j;—r) >> a^;, which defines the internal 
region. Likewise [1 — -F(r)] — >■ 1 for (r — R^) » a^, which defines the external region. 
Thus Rx forms an approximate border between the internal and the external regions, and 
if Rx is close to the nuclear surface then the two regions will correspond to the high density 
and low density regions of nuclear matter, respectively. The fixed 6o in the tp potential was 
replaced by 

6o ^ 5o F(r) + 6o [1 - F{r)] (48) 
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FIG. 22: Contributions to the of K atomic shifts (left) and widths (right) for the deep density- 
dependent potential F from Ref. [108] and for the shallow chirally-based potential from Ref. 132]. 



where the parameter represents the low-density interaction and the parameter Bq repre- 
sents the interaction inside the nucleus. This division into regions of high and low densities 
is meaningful provided is close to the radius of the nucleus and is of th e order of 
0.5 fm. This is indeed the found in global fits to kaonic atom data |108[ |. We note 



that unlike with pionic and antiprotonic atoms, the dependence of kaonic atom fits on the 
rms radius of the neutron distribution is marginal. 

Figure [21] (left) shows, as an example, the real part of the best-fit potential for ^^Ni 
obtained with the various models discussed above, i.e. the simple tp model and its DD 
extension, and the geometrical model F, with the corresponding values of fo^' 65 data 



points in parentheses. Also shown, with an error band, is a Fourier-Bessel (FB) fit |10l| that 
is discussed below. We note that, although the two density-dependent potentials marked 
DD and F have very different parameterizations, the resulting potentials are quite similar. 
In particular, the shape of potential F departs appreciably from p(r) for p(r)/po < 0.2, 
where the physics of the A (1405) is expected to play a role. The density dependence of the 
potential F provides by far the best fit ever reported for any global i^'~-atom data fit, and 
the lowest value as reached by the model-independent FB method. On the right-hand 
side of the figure are shown the individual contributions to of the shifts for the deep F 
potential and the shallow tp potential. Figure [22] shows comparisons between values of 
the shifts and of th e widths for the F potential and the yet shallower chirally-based potential 
of Baca et al. 



132 



It is self evident that the agreement between calculation and experiment 
is substantially better for the deep F potential than for the shallow chiral potentials. 

The questio n of how well the real part of the iir~-nucleus potential is determined was dis- 



cussed in Ref. |101| . Estimating the uncertainties of hadron-nucleus potentials as function of 
position is not a simple task. For example, in the 'tp' approach the shape of the potential is 
determined by the nuclear density distribution and the uncertainty in the strength parame- 
ter, as obtained from x"^ fits to the data, implies a fixed relative uncertainty at all radii, which 
is, of course, unfounded. Details vary when more elaborate forms such as 'DD' or 'F' are 
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FIG. 23: Functional derivatives of kaonic atoms with respect to the fully complex (Comp, 
dashed lines) and real (Re, solid lines) potential as function of r/, where r = Rc + rjac, with R, 
and Oc the radius and diffuseness parameters, respectively, of th e cha rge distribution. Results are 
shown for the tp potential and for the t{p) 'F' potential of Ref. jlOg ] obtained from global fits to 
kaonic atom data. 



used, but one is left essentially with analytical continuation into the nuclear interior of poten- 
tials that might be well-determined only close to the nuclear surface. 'Model-independent' 
methods have been used in analyses of elastic scattering data for various projectiles to 
alleviate this problem. However, applying e.g. the Fourier-Bessel (FB) method in global 
analyses of kaonic atom data, one ends up with too few terms in the series, thus making the 
uncertainties unrealistic in their dependence on position. This is illustrated in Fig. [5T]by the 
'FB' curve, obtained by adding a Fourier-Bessel series to a 'tp' potential. Only three terms 
in the series are needed to achieve a of 84 and the potential becomes deep, in agreement 



with the other two 'deep' solutions. The error band obtained from the FB method is, 
nevertheless, unrealistic because only three FB terms are used. However, an increase in the 
number of terms is found to be unjustified numerically. 

The functional derivative (FD) method for identifying the radial regions to which exotic 
atom data are sensitive is described in detail in Sec. IHID 4[ This method was applied in 



Ref. [101| to the F and tp kaonic atom potentials and results are shown in Fig. [23] where rj 



is a global parameter defined by r = i?c + V^^c, with and Oc the radius and diffuseness 
parameters, respectively, of the charge distribution. From the figure it can be inferred that 
the sensitive region for the real tp potential is between t] = —1.5 and rj = 6 whereas for 
the F potential it is between t] = —3.5 and rj = 4. Recall that rj = —2.2 corresponds 
to 90% of the central charge density and t] = 2.2 corresponds to 10% of that density. It 
therefore becomes clear that within the tp potential there is no sensitivity to the interior of 
the nucleus whereas with the t{p) 'F' potential, which yields greatly improved fit to the data, 
there is sensitivity to regions within the full nuclear density. The different sensitivities result 
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FIG. 24: Calculated energies of K~ atomic states in ^'^^Pb. The lowest energy for each I value 
corresponds to n = / + 1. The bars represent the widths of the states. 



from the potentials themselves: for the tp potential the interior of the nucleus is masked 
essentially by the strength of the imaginary potential. In contrast, for the F potential not 
only is its imaginary part significantly smaller than the imaginary part of the tp potential 



108| but also the additional attraction provided by the deeper potential enhances the atomic 



wavef unctions within the nucleus [1] thus creating the sensitivity at smaller radii. As seen 
in the figure, the functional derivative for the complex F potential is well approximated by 
that for its real part. 

It is concluded that optical potentials derived from the observed strong-interaction effects 
in kaonic atoms are sufficiently deep to support strongly-bound antikaon states, but it does 
not necessarily imply that such states are sufficiently narrow to be resolved unambiguously 
from experimental spectra^ Moreover, the discrepancy between the very shallow chirally 



motivated potentials [2J, [llOj, the intermediate potentials of depth around 100 MeV [26| 
and the deep phenomenological potentials of type 'F' remains an open problem. It should 
also be kept in mind that these depths relate to K potentials at threshold, whereas the 
information required for ^-nuclear quasibound states is at energies of order 100 MeV below 
threshold. 



E. Deeply bound K atomic states 

Somewhat paradoxically, due to the strong absorptive imaginary part of the /^"-nucleus 
potential, relatively narrow deeply-bound atomic states are expected to exist which are quite 
independent of the real potential. Such states are indeed found in numerical calculations as 
can be seen in Fig. [2l] where calculated binding energies and widths of atomic states of K~ 
in 208p]^ g^j,g shown for several /-values, down to states which are inaccessible via the X-ray 
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FIG. 25: Saturation of width F for the 2p 'deeply bound' K atomic state in ^'^^Pb as function of 
absorptivity, Im^o, for Re6o = 0.62 fm (left) and wavefunctions for this state (right). 



cascade. For ^°^Pb, the last observed K~ atomic circular state is the 7i, corresponding to 
/ = 6. The general physics behind this phenomenon is similar to that responsible for the 
deeply-bound pionic atom states, although there are differences in the mechanism. For a 
Schrodinger equation the width of a state is given exactly by Eq. (l23l) and if a normalized 
atomic wavefunction is expelled from the nucleus, then small widths are expected due to 
the reduced overlap between the atomic wavefunction and the imaginary po tential. ( For the 
KG equation there are small changes in the expression Eq. fl23|) . see Refs. jl33l . 134| -) The 
mechanism behind the pionic atom deeply bound states is simply the repulsive real part of 
the s-wave potential. In contrast, phenomenological kaonic atom potentials are attractive, 
but the strengths of the imaginary part of the potential are such that the decay of the 
wavefunction as it enters the nucleus is equivalent to repulsion, resulting in narrow atomic 
states due to the reduced overlap as discussed above. It is seen from Fig. [2l] that there 
is a saturation phenomenon where widths hardly increase for / < 2, contrary to intuitive 
expectation. The repulsive effect of sufficiently strong absorption is responsible for the 
general property of saturation of widths of atomic states and saturation of reaction cross 
sections above threshold, observed experimentally for antiprotons (see Sec. IVIFI ) 

The left-hand side of Fig. [25] shows the saturation of widths as function of Im6o for the 
2p state of kaonic atoms of ^°^Pb. For small values of Imfeg the calculated width increases 
linearly but already at 20% of the best-fit value of 0.9 fm saturation sets in and eventually 
the width goes down with further increase of the absorption. Note that the real part of the 
binding energy, represented here by the strong interaction level shift e, is hardly changing 
with Im6o- The right-hand side of Fig. [25] shows radial wavefunctions for the 2p atomic 
K~ state in ^°^Pb for several combinations of potentials. The dashed curve marked 'Coul' 
is for the Coulomb potential only, and with a half-density radius for ^''^Pb of 6.7 fm it 
clearly overlaps strongly with the nucleus. Adding the full complex optical potential the 
solid curve marked 'Comp' shows that the wavefunction is expelled from the nucleus, and 
the dotted curve marked 'Im' shows that this repulsion is effected by the imaginary part 
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of the potential. Clearly the overlap of the wavefunction with the nucleus is dramatically 
reduced compared to the Coulomb-only situation. An interesting phenomenon is displayed 
by the dot-dashed curve marked 'Re'. It shows the wavefunction when the real potential 
is added to the Coulomb potential, demonstrating significant repulsion of the wavefunction 
by the added attractive potential. The explanation for this bizarre result is provided by 
the three small peaks inside the nucleus which are due to the orthogonality of the atomic 
wavefunction and strongly-bound K~ nuclear wavefunctions having the same /-values. This 
extra structure of the wavefunction in the interior effectively disappears when the imaginary 
potential is included. 



F. Deeply bound K nuclear states in light nuclei 



No saturation mechanism holds for the width of ^-nuclear states which retain very good 
overlap with the potential. Hence, the questions to ask are (i) whether it is possible at all 
to bind strongly K mesons in nuclei, and (ii) are such quasibound states sufficiently nar- 
row to allo w obs ervation and identification? The first question was answered affirmatively 
by Nogami |l35 | as early as 1963 arguing that the K~pp system could acquire about 10 
MeV binding in its / = 1/2, L = S = state. Yamazaki and Akaishi, using a complex 
energy-independent KN potential within a single-channel K~pp calculation [136i] , reported 
a binding energy i? ~ 50 MeV and width F ~ 60 MeV. Preliminary results of revised Anti- 
symmetrized Molecular Dynamics (AMD) calculations by Dote and Weise jl24 . 137l | which 
implicitly account for KN — vrE coupling agree on -B but estim ate F ~ 100 MeV. Coupled- 
channel KNN — ir'EN Faddeev calculations [138 . 139l . 140, 141 | of K~pp have confirmed this 
order of magnitud e of bindi ng, B ~ 55-75 MeV, differing on the width; the calculations by 
Shevchenko et al. 138 . [139I give large values, F ~ 100 MeV, for the mesonic width. These 
Faddeev calculations overlook the KNN YN coupling to nonmesonic channels which 
are estimated to add, conservatively, 20 MeV to the overall width. Altogether, the widths 
calculated for the K~pp quasibound state are likely to be so large as to make it difficult to 



identify it experimentally [142 



The current experimental and theoretical i ntere st in i^-nuclear bound states was triggered 
back in 1999 by the suggestion of Kishimoto 1431 ] to look for such states in the nuclear re- 



145[ | who suggested to look 
ttS decay 



action {K~,p) in fiight, and by Akaishi and Yamazaki [144, 
for a KNNN 1 = state bound by over 100 MeV for which the main KN 
channel would be kinematically closed. In fact , Wycech had conjectured that the width of 
such states could be as small as 20 MeV [l46 |. Some controversial evidence for relatively 
narrow states was presented initially in (A'stopiP) (-^stop;'"') reac tions on ^He (KEK-PS 

^-nuclear states 



E471) [33, [381 but has recently been withdrawn (KEK-PS E549/570) [147 , 
were also invoked to explain few statistically-weak irregularities in the neutron spectrum of 
the {K~, n) in-fiight reaction on ^^O (BNL-AGS, parasite E930) j39|, but subsequent {K~ ,n] 



and {K-,p) reactions on ^^C at pi^h = 1 GeV/c (KEK-PS E548) |l48[ | have not disclosed 



any peaks beyond the appreciable strength observed below the ^-nucleus threshold. Ongo- 
ing experiments by the FINUDA spectrometer collaboration at DA$NE, Frascati, already 
claimed evidence for a relatively broad K^pp deeply bound state {B ~ 115 MeV) in K^p 
reactions on Li and ^^C, by observing back-to-back Ap pairs from the decay K~pp — ^ Ap (40|, 
but these pairs could naturally arise from conventional absorption processes at rest when 
final-state interaction is taken into account 
^Li has also been recently observed 



149j. Indeed, the K^^^ppn — S p reaction on 



150|. Another recent claim of a very narrow and deep 



45 



K pp state {B ~ 160 MeV, F ~ 30 MeV) is also based on observing decay Ap pairs , using 
p annihilation data on ^He from the OBELIX spectrometer at LEAR, CERN jl51 |. The 



large value of B^-pp over 100 MeV conjectured by these experiments is at odds with all 
the few-body calculations of the K~pp system listed above. One cannot rule out that the 
pairs assigned in the above analyses to K~pp decay in fact result from nonmesonic de- 
cays of different clusters, say the KNNN 1 = quasibound state. A definitive study of 
the K~pp quasibound state (or more generally {-ft'[A^A^]/=i}/=i/2) could be reached through 
fully exclusive formation reactions, such as: 

ir-+3He ^ n + {K[NN]i=^]j=y2,i.=+i/2, p+{K[NN]i=^]j=y2,i.=_y2, (49) 

the first of which is scheduled for day-one experiment in J- PARC jl52l | . Finally, preliminary 
evidence for a KNNN 1 = state with 5 = 58±6 MeV, F = 37±14 MeV has been recently 
presented by the FINUDA collaboration on ®Li by observing back-to-back Ad pairs [l53']. It 
is clear that the issue of K nuclear states is far yet from being experimentally resolved and 
more dedicated, systematic searches are necessary. 



G. RMF dynamical calculations of K quasibound nuclear states 



In this model, spelled out in Refs. [l08l . Il54j . Il55| . the (anti)kaon interaction with the 



nuclear medium is incorporated by adding to Cjy the Lagrangian density Ck- 

Ck = V*^KV^'K - mlKK - g,KmK<yKK . (50) 

The covariant derivative T>^ = + igujK^n describes the coupling of the (anti)kaon to 
the vector meson u. The(anti)kaon coupling to the isovector p meson was found to have 
negligible effects. The K meson induces additional source terms in the equations of motion 
for the meson fields a and Uq. It thus affects the scalar S = q^nO' and the vector V = Qujn^o 
potentials which enter the Dirac equation for nucleons, and this leads to rearrangement or 
polarization of the nuclear core, as shown on the left-hand side of Fig. [26] for the calculated 
average nuclear density P = \ ! P^dr as a function of B^- for K~ nuclear Is states across 
the periodic table. It is seen that in the light K~ nuclei, p increases substantially with 
Bx- to values about 50% higher than without the K. The increase of the central nuclear 
densities is bigger, up to 100%, and is nonnegligible even in the heavier K~ nuclei where 
it is confined to a small region of order 1 fm. Furthermore, in the Klein-Gordon equation 
satisfied by the K, the scalar S = gaK<7 and the vector V = —Quk^q potentials become 
state dependent through the dynamical density dependence of the mean-field potentials S 
and as expected i n a R MF calculation. An imaginary ImVjf ~ tp was added, fitted to 
the K~ atomic data jl07l |. It was then suppressed by an energy-dependent factor f{Bj^), 
considering the reduced phase-space for the initial decaying state and assuming two-body 
final-state kinematics for the decay products in the KN tcY mesonic modes (80%) and 
in the KNN YN nonmesonic modes (20%). 

The RMF coupled equations were solved self-consistently. For a rough idea, whereas 
the static calculation gave = 132 MeV for the K~ Is state in ^^C, using the 

values g^j^i q'^k^ corresponding to the i^~-atom fit, the dynamical calculation gave 
B]^^ = 172 MeV. In order to scan a range of values for -B^-, the coupling constants 
QaK and Qujk were varied in given intervals of physical interest. 

Beginning approximately with ^^C, the following conclusions may be drawn: 
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FIG. 26: Dynamically calculated average nuclear density p (left) and widths F^ - (r ight) of Is 
i^^-nuclear states in the nuclei denoted, as function of the Is K~ binding energy 108]. 



For given values of gaK-, QujK^ the K binding energy Bj^ saturates as function of A, 
except for a small increase due to the Coulomb energy (for K~). 

The difference between the binding energies calculated dynamically and statically, 
B'^^ — B^^^^, is substantial in light nuclei, increasing with Bj^ for a given value of A, 
and decreasing monotonically with A for a given value of Bj^. It may be neglected only 
for very heavy nuclei. The same holds for the nuclear rearrangement energy B^^' — Bj^ 
which is a fraction of B'^^ — B^^^. 

The width F^(i?^) decreases monotonically with A, according to the right-hand side 
of Fig. [26] which shows calculated widths F;^- as function of the binding energy B^- 
for Is states in ^-^C, ^l°Ca and ^^-''^Pb. The dotted line shows the static 'nuclear- 
matter' limit corresponding to the K~-a,tom. fitted value Im6o = 0.62 fm and for 
p{r) = po = 0.16 fm~^, using the same phase-space suppression factor as in the 
'dynamical' calculations. It is clearly seen that the functional dependence F/^-(i?x-) 
follows the shape of the dotted line. This dependence is due primarily to the binding- 
energy dependence of the suppression factor / {Bk-) which falls off rapidly until B^- ~ 
100 MeV, where the dominant KN vrS gets switched off, and then stays rather flat 
in the range B^- ~ 100-200 MeV where the width is dominated by the KNN YN 
absorption modes. The widths calculated in this range are considerably larger than 
given by the dotted line (except for Pb in the range B^- ~ 100-150 MeV) due to the 
dynamical nature of the RMF calculation, whereby the nuclear density is increased by 
the polarization effect of the K~ . Adding perturbatively the residual width neglected 
in this calculation, partly due to the KN vrA secondary mesonic decay channel, a 
representative value for a lower limit T ^ > 50 ± 10 MeV holds in the binding energy 
range B^- ~ 100-200 MeV. Fig. [27] shows the effect of splitting the 80% mesonic 
decay width, previously assigned all to vrS absorption channels, between ttE (70%) 
and ttA (10%), and also of simulating the 20% nonmesonic absorption channels by 
a p^ dependen ce c ompared to ImV^ ~ tp used by Mares et al. [108]. These added 
contributions |l55l | make the lower limit F^ > 50± 10 MeV a rather conservative one. 
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nuclear state in j^'^Ca for vrS + ttA com- 



FIG. 27: Dynamically calculated widths of the Is K 
pared to vrS final mesonic absorption channels, and for compared to p dependence of t he fi nal 
nonmesonic absorption channels, as function of the K~ binding energy, from Gazda et al. 155 1. 



V. K+ MESONS 



Overview of the X^-nucIeus interaction 



There are obviously no exotic atoms to provide information on the i^"'"-nucleus inter- 
action below threshold. Nevertheless, it has been found that for pions and antiprotons the 
optical potentials cross smoothly from the atomic into the scattering regime, and therefore 
studies of the /^"'"-nucleus interaction at low energies above threshold are relevant to the 
general topic of medium modifications of the interaction. In fact, mesons provide a 
clear example of such modifications. 

The K^N interaction below the pion-production threshold is fairly weak and featureless 
and this merit has motivated past suggestions to probe nuclear in-medium effects by st udying 



scattering and reaction processes with beams below 800 MeV/c; see Ref. |l56j | for an 
early review. The insufficiency of the impulse-approximation tfreeP form of the - nucleus 
optical potential, where tfree is the free-space K^N t ma trix, was somewhat of a surprise 
already in the 1980s. Limited total cross-section data 157l | on carbon, and elastic and 
inelastic differential cross section data [158| on carbon and calcium showed problems with 

the tp potential, particularly with respect to its reaction content ('reactivity' below); In 

order to account for the increased reactivity in - nucleus interactions, Siegel et al. jl59l 



and later Peterson [160[ | suggested that nucleons 'swell' in the nuclea r me dium, primarily 



by increasing the dominant hard-core Su phase shift. Brown et al. [161| suggested that 



the extra reactivity was due to the reduced in-medium mass es of exchanged vector mesons, 
and this was subsequently worked out in detail in Ref. 
reactivity in - nucle us intera ctions was discussed 
exchange- current effects jl63 . 164 . 



162|. Another source for increased 
the 1990s and is due to meson 
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Some further experimental progress was made during the early 1990s, consisting mostly 
of measuring attenuation cross sections in transmission experiments at the BNL-AGS on 
deuterium and several other nuclear targets in the momentum range piab = 450-740 MeV/c 
165 . 1661 . 167 | and extracting i^'"^-nucleus total cross sections. The same transmission data 



were then reanalyzed to extract in addition total reaction cross sections [41], which are less 
dependent on the potential used in converting measured attenuation into cross sections and 
eventually self-consistent final values of integral cross sections (reaction and total) on 
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^Li, ^^C, ^'^Si and ^°Ca were pubhshed in Refs. [4 

gave clear evidence for the density dependence of the increased reactivity suggested by the 

a. 



These integral cross-section data 
eased reactivity suggested by the 
earlier data, as demonstrated in Fig. [2H] from Ref . [42] . Plotted in the figure are ratios of 
experimental to calculated integral cross sections for ^^C, ^^Si and ^°Ca, where the calculated 
cross sections use a tp potential fitted to the ^Li data. These ratios, for the denser nuclei, 
deviate considerably from the value of one in a way which is largely independent of the beam 
momentum. 

Other measure ments during the 1990s include quasielastic scattering on several tar- 



gets at 705 MeV/c [168| and new meas urern ents of elastic and inelastic differential cross 



sections on ^Li and at 715 MeV/c ^M, 



These data and 



further analyzed in Ref. 
analyses lent support to the substantial medium modifications demonstrated above on the 
basis of studying integral cross sections. By the late 1990s, experimentation in - nuclear 
physics had subsided, and with it died out also theoretical interest although the problems 
with medium modifications of the interaction remained, as is shown be low. Theoretical in - 
terest in - nuclear physics to some extent has been revived recently |l71 . 172, 173 . 174], 
particularly in connection with possible contributions due to G"*" pentaquark degrees of free- 
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dom, as is also discussed below. 



B. Kaon-nucleus optical potential 

The starting form adopted for the kaon-nucleus optical potential Vopt, following Eq. (jl]) 
of Sec. II Bt is the simplest possible tp form: 

24^iKpt(r) = -47rF^&0P(r) , (51) 

where e^^\ is the center-of-mass (cm.) reduced energy, 

= E;' + E'/ (52) 

in terms of the cm. total energies for the projectile and target respectively, and 



M{Ea + Ep 



(53) 



is a kinematical factor resulting from the transformation of amplitudes between the KN and 
the - nucleus cm. systems, with M the free nucleon mass. Ma the mass of the target 
nucleus and ^/s the total projectile-nucleon energy in their cm. system. The parameter 
bo in Eq. flSTl) reduces in the impulse approximation to the (complex) isospin-averaged KN 
scattering amplitude in the forward direction. For ^Li and for ^^C the modified harmonic 
oscillator (MHO) form was used for the nuclear densities whereas for ^^Si and for ^°Ca the 
two-parameter Fermi (2pF) form was used and minor changes were made to parameters of 
the neutron density to check sensitivities to p„. The data base for the analysis were the 32 
integral cross sections for on ^Li, ^^C, ^'^Si an d ^°Ca fro m Ref. (43| . 

Fits to the integral cross sections were made 171 . 172l | separately at each of the four 



momenta, varying the complex parameter bo, and the results are summarized in Table IVnl 
marked as tp for each momentum. From the values of P^r point it is seen that the 
fits are unacceptably poor and the resulting Re6o and Im6o disagree with the correspond- 
ing fre e K '^N values (marked as tfreeP and derived from the KN phase shifts as given by 
SAID jl75l |). The discrepancies are particularly noticeable for Im6o, which are determined to 
good accuracy. Evidently the experimental results indicate significant increase in reactivity, 
as mentioned in Sec. IV A[ 

The obvious next step is to effectively make bo density dependent by introducing higher 
powers of the density, such as 

bo p(r) ^ bo p(r) + B p^{r) , (54) 
where bo th parameters bo and B are to be determined from fits to the data. The results 



17ll . Il72[ | are also shown in Table IVIIl marked as Eq. flM|) , where Refeo was kept fixed at 
its respective free KN value. The improvement in the fits to the data is evident from the 
reduction of the values of x^? but except for the lowest incoming momentum the quality 



of the fits suggests that something is still missing. Guided by earlier analyses |42l . |43| that 
achieved much improved fits by introducing the average nuclear density p for each target 
nuclide, 



\ I P'dr , (55) 
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TABLE VII: Fits to the eight - nuclear integral cross sections [43| at each of the four laboratory 
momenta piab (in MeV/c), using different potentials. 



Plah 


Vovt 


Refto (fill) 


Im6o(fm) 


ReS(fm'^) 


ImS(fm4) 




488 


tp 


-0.203(26) 


0.172(7) 






16.3 






-0.178 


0.153 










Eq.dMD 


-0.178 


0.122(5) 


0.52(20) 


0.88(8) 


1.18 




Eq (l56l) 


-0.178 


0.129(4) 


0.17(11) 


0.62(6) 


0.27 


531 


tp 


-0.196(39) 


0.202(9) 






56.3 






-0.172 


0.170 










Eq.dMI) 


-0.172 


0.155(14) 


1.79(46) 


0.72(27) 


7.01 




Eq (1561) 


-0.172 


0.146(5) 


0.46(21) 


0.78(7) 


3.94 


656 


tp 


-0.220(50) 


0.262(12) 






54.9 






-0.165 


0.213 










Eq.dMI) 


-0.165 


0.203(18) 


1.66(80) 


0.89(36) 


7.24 




Eq.dMI) 


-0.165 


0.204(5) 


2.07(19) 


0.77(7) 


0.32 


714 


tp 


-0.242(53) 


0.285(15) 






67.7 






-0.161 


0.228 










Eq.dMI) 


-0.161 


0.218(24) 


1.40(95) 


1.10(48) 


9.3 




Eq.dMD 


-0.161 


0.218(6) 


1.51(43) 


0.97(9) 


1.24 



Eq. fl54|) was replaced by the following ansatz 171 , 172l 



bo p{r) bo p{r) + B p p{r) . (56) 

The results of this prescription are also shown in Table IVIIl marked as Eq. fl56l) and it is 
clear that the fits to the data are very good. This use of the average nuclear density singles 
out the ^Li target from the other three targets, due to its average density p being close to 
50% of the corresponding values for the other targets. This is a purely phenomenological 
observation without (as yet) any theory behind it. 

In order to further test the picture that emerges from the analysis of th e integral cross 
sections for the - nucleus interaction, the analysis was repeated jl72l | including also 



differential cross sections for the elastic scattering of by some of the target nuclei. Fits 
were made to the combined integral and differential cross sections at 714 Me V/c consisting 
of the eight integral cross sections and the 17 differential cross sections from Ref. usmg 



the potentials of either Eq. (15^ or Eq. (!56|) . Figure [291 shows that again the pp version of 
the potential Eq. (l56l) is preferred and the fits to the differential cross sections are good. 
The potential parameters obtained from the fits to the combined integral and differential 
cross sections agree, within uncertainties, with the corresponding values in Table IVIIl 

Prior to discussing in the next section the reactive content of the above forms of density- 
dependent - nucleus optical potentials, it is worth noting that the splitting of ImV^pt 
in Table IVIIl into its two reactive components Imbo and ImB appears well determined by 
the data at all energies, and perhaps is even model independent, particularly for the pp 
version Eq. fl56l) of the optical potential for which very accurate values of \mbo are derived. 
These values of Im6o are close to, but somewhat below the corresponding free-space values, 
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cm. angle (deg.) cm. angle (deg.) 

FIG. 29: Comparison between m easu red differential cross sections for elastic scattering at 
Piab = 715 MeV/c on ^Li and [l70( and best-fit calculations using Eq. ([5^ (dashed lines) and 
Eq. ([55]) (solid lines). 



a feature which is observed in calculations which replace the tp form hy gp where nuclear 



phase space effects are considered explicitly [173|. The values derived for ImB are roughly 



independent of the form of the piece added to tp, AKpt, whether Eq. fl54l) or Eq. fl56l) are 
used to derive these values from the data. In contrast, the two components of ReKpt are 
correlated strongly when Re6o is also varied, largely cancelling each other into a resulting 
poorly determined ReV^pt. For this reason, it will appear difficult to offer any conclusive 
model for the physics underlying the real part of Al^pt. 



C. absorption cross sections 

The effect of AVopt, within the improved fits to the i^'^-nuclear integral cross sections, 
is demonstrated in the upper part of Fig. [30] for the reaction cross sections per nucleon 
aii/A at 488 MeV/c, where the calculated cross sections using a best-fit to optical potential 
(dashed line) are compared with the experimental values listed in Ref. [43|. The tp fit fails 
to reconcile the ^Li data with the data on the other, denser nuclei. If ^Li is removed from 
the data base, then it becomes possible to fit reasonably well the data for the rest of the 
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FIG. 30: Data and calculations 17l[ | for reaction cross sections per nucleon (aR/A) at piab = 
488 MeV/c are shown in the upper part. Calculated absorption cross sections per nucleon 
(cabsM) are shown in the lower part, see text. 



nuclei, but the rise in Im^o with respect to its free-space value is then substantially higher 
than that for the tp potential when ^Li is included. At the higher energies, tp fits which 
exclude ®Li are less successful than at 488 MeV/c, while also requiring a substantial rise in 
Imfeo) which means increased values of the in-medium KN total cross sections with respect 
to the corresponding free-sp ace v alues. This has been observed also in a - nucleus 
quasifre e-sca ttering analysis |l76 |. for mesons incident on C, Ca, Pb at piab = 705 



MeV/c |168| . Also shown in the upper part of Fig. [30], by the solid line marked tp + AVopt, 
are calculated reactio n cro ss sections at 488 MeV/c using Eq. (!56|) . This is a very good fit. 
Recently, Tolos et al. 173l | have demonstrated that a similarly substantial improvement in 
the reproduction of reaction cross sections could be achieved microscopically by coupling 
in pentaquark degrees of freedom. It is tempting to assume that the effects of absorbing 
mesons into a pentaquark configuration are given by the difference between the solid 
curve and the dashed curve in Fig. [20] for nuclear targets heavier than ^Li. However, for a 
quantitative estimate of the absorption cross sections one needs to do a more explicit 
calculation. 

In close analogy to analyses of pionic atoms and of low-energy pion-nuclear scattering 
reactions in which the parameter Bq (cf. Eq. f[T5l) ) is related to vr" nuclear absorption 
processes on two and on more nucleons, the additional piece Al^pt due to the nonzero value 
of the parameter B in Eqs. (I54|56p represents nuclear absorption into B"*" - nuclear 
final states. Here 0"*", with mass Mq+ ^ 1540 MeV, is an hypothetical S = +1 'exotic' 
pentaquark bar yon, searches for which have not provided conclusive evidence (for a recent 



review see Ref . 177| ) . The abnormally small upper limit < 1 MeV for the width of the 
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6+ pentaquark deduced from some of these searches imphes a neghgible 9+ KN couphng, 
but this does no t hrn it the couphng Q^N — » KNN which is related to virtual decays such 
as G"*" — i> ttKN jl78 |. Gal and Friedman jl7ll . 172] estimated the nuclear absorption cross 



sections of mesons by using two slightly different versions of the distorted-wave Born 
approximation: 

^^""'^ -'I;/ Ini(AKpt(r)) |^S:l„,,=o)(r)r dr , (57) 



a. 



abs 



and 



a 



{K+) 
abs 



- ^1 Im(AV;pt(r)) |^(+)(r)p dv , (58) 

where the distorted waves '^'^(^Vo t=o) calculated discarding AKpt. 

Calculated absorption cross sections per target nucleon at piab = 488 MeV/c are shown 
in the lower part of Fig. [30] for the fit using Eq. (|56l) for Vo^t in Table I VIII The triangles are 
for expression fl57|) and the solid circles are for expression fl58l) . The error bars plotted are 
due to the uncertainty in the parameter Imi?. It is seen that these calculated absorption 
cross sections, for the relatively dense targets of ^^C, ^^Si and ^''Ca, are proportional to the 
mass number A, and the cross section per target nucleon due to \m.B 7^ is estimated as 
close to 3.5 mb. Although the less successful Eq. (IMl) gives cross sections larger by 40% 
at this particular incident momentum, this value should be regarded an upper limit, since 



the best-fit density-dependent potentials of Refs. [42|, |43| yield values smaller th an 3 .5 mb 
by a similar amount. The experience gained from studying vr- nuclear absorption |l79 | leads 

to the conclusion that aahs{K~^ N N) is smaller than the extrapolation of cr^^ ^ /A in Fig. [30] 
to A = 1, and since the KN interaction is weaker than the ttN interaction one expects a 
reduction of roughly 50%, so that aa_^,s{K~^NN) ~ 1-2 mb. 

In Fig. [301 the considerably smaller absorption cross section per nucleon calculated for the 
relatively low-density ®Li nucleus suggests a cross section of order fraction of millibarn, in a 
possible missing-mass search based on observing the final proton in the two-body reaction 
K~^d — ^ Q~^p. This cross section is not expected to exhibit marked resonance behavior near 
Piah ~ 440 MeV/c, which corresponds to the G"'"(1540) resonance assumed rest mass, even 
if B"*" is very narrow. For the heavier nuclear targets too, the assignment of the excess 
reactivity observed in i^^-nuclear cross sections as due to = -|-1 pentaquark degrees of 
freedom does not require the existence of a narrow KN resonance. It only assumes that 
pentaquark degrees of freedom are spread over this energy range with sufficient spectral 
strength. For nuclear targets other than deuterium, given the magnitude of the K~^ nuclear 
absorption cross sections as reviewed here, {K~^,p) experiments could prove useful. This 
reaction which has a 'magic momentum' about pi^h ~ 600 MeV/c, where the 0"'"(154O) is 
produced at rest, is particularly suited to study bound or continuum states in h yponuclei 
{S = +1 nuclei according to the terminology suggested by Alfred Goldhaber jl80l |). 

In conclusion, the available K~^ nuclear cross section data at puh ~ 450-800 MeV/c 
reveal substantial reactivity beyond that produced by the impulse approximation, or for that 
purpose by any effective tp form of the K~^ optical potential. It was shown that this extra 
reactivity may be explained by adding a two-nucleon absorption channel K^nN —>■ Q~^N, 
where the degrees of freedom need not materialize within a narrow energy bin. This 
provides a density-dependent mechanism that couples in S" = +1 pentaquark degrees of 
freedom in a way which is insensitive to the width of their spectral distribution ^ITll . While 
there is no firm support at present for this conjecture from other phenomenological sources, 
a robust experimental program of measuring low-energy K~^d and K^ - nuclear scattering 
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and reaction cross sections in the range piab ~ 300-800 MeV/c, and particularly about 400 
Mev/c, would be extremely useful to decide whether or not S = +1 pentaquark degrees of 
freedom are involved in - nuclear dynamics. 



VI. ANTIPROTONS 



A. Overview of the p-nucleus potential 

In line with the other types of exotic atoms, the interaction of antiprotons with nuclei 
at threshold is described in terms of an optical potential, which in the simplest tp form is 
given by 

2/xKpt(r) = -47r(l + A^i_l)[5p(p^ + p^) + 5^(p„ _ p^)] , (59) 

where /i is the reduced mass of the p, pn and pp are the neutron and proton density dis- 
tributions normalized to the number of neutrons and number of protons Z, respectively, 
A = N + Z, and M is the mass of the nucleon. The factor {A — I) /A above, which was 
omitted from the potential for pions, is included here due to the larger mass of the p. Be- 
cause of the large cross section for annihilation of p on a single nucleon, the interaction is 
expected to be dominated by the imaginary part of the potential and the absorption of p 
is expected to take place at the extreme surface regions of the nucleus. As a result it is 
unlikely that p atoms will provide information on the potential deep into the nucleus and 
the above simplest tp form is a useful starting point for analyzing antiprotonic at om d ata. 



Previous attempts to add to the potential a p-wave term or non linear terms [H, II8I 



are 



not found to be justified in a phenomenological approach when the overall picture is con- 
sidered, respecting also constraints satisfied by neutron density distributions, as described 
below. More specifically, we f i nd that a n imaginary part of a p-wave potential compatible 
with the Paris potential |l82 . 183l . 184 1 could be accommodated, but then Imfeg is found 



to be incompatible with the Paris potential. These remarks apply particularly to analyses 
based on the high-quality data of the PS209 collaboration, which is the basis for the present 
analysis. However, an isovector term bi {pn — pp) is included in Eq. f l59|) because the present 
data base is rich with groups of isotopes of the same element. 

Proton densities for the above potential are taken, as before, from the known charge 
distributions |10| by unfolding the finite size of the proton. For the neutrons it is, again, a 
matter of choosing an adequate model, that will be in line with the bulk of information on 



neutron densities [16[ . The importance of the shapes of neutron dens ity d istributions pn was 



realized long ago, when single-particle densities (SP) were found [H. Il81 | to produce better 



fits to the data of p atoms compared to fits based on the two-parameter Fermi (2pF) shape, 
because of the sensitivity of p atom data to the extreme outer reaches of the nucleus. An 
alternative to SP densities is to use simple parameterizations such as the 2pF form for p„ 
and to accommodate different shapes, relative to the protons, by taking the 'skin' or the 
'halo' version for p„ fl6|], or their average. In what follows we adopt the latter approach for 
global fits to p data because in any case the SP densities are not particularly suitable for 
nuclei far removed from closed shells and because we are interested in average properties. 
We therefore use the approach of Sec. IIII D 21 with the parameterization of r„ — given by 
Eq. ©. 
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B. Antiprotonic atom data 



As mentioned in Sec. HIl the experimental situation with p atoms has changed significantly 
in the last decade with the publication by the PS209 collaboration [48| of high-quality data 
for several sequences of isotopes along the periodic table. This set of X-ray data has greater 
accuracies compared to older data used in earlier analyses, and thanks to the full coverage of 
the periodic table by these new data we do not mix in the present analysis old data with the 



new results but use only th e PS209 r esults |48| . including the revised experimental results 



for Cd, Sn and Te isotopes |l85 . [l86l obtained after correcting for E2 resonance effects. In 
what follows we address only spin- averaged quantities for antiprotonic atoms within a global 
approach to the hadron-nucleus interaction. 

Strong interaction effects in antiprotonic atoms are reported j48| as level shifts and widths 
for the lowest levels reached in the X-ray cascade and as 'upper' level widths deduced from 
the yields of the transitions, based on intensity ratios and calculations of the atomic cascade 



process. These yield values had been converted in Ref. [48|| into upper level widths with 
the help of the calculated radiation widths. However, it is easy to see that values for the 
deduced upper level widths may be different from the corresponding values calculated for 
the yields. As the yields are the experimentally determined quantities, we have converted 
the quoted upper level widths back to transition yields, and used these in the global fits. 
The X-ray data used in the present analysis are for the following nuclear targets: ^^'^^O, 

40,42,43,44,48Qg^ 54,56,57,58pg 58,60,62,64j^j 90,962^. 106,116q^ 112,116,120,124g^ 122, 124,126, 128,130rpg 

208pb, a total of 90 data points Hi . 



In addition to the conventional method of studying strong interaction effects by observa- 
tion of X-ray emission from exotic atoms, there is for antiprotons a radi oche mical method 
which is capable of providing information on the absorption of p by nuclei jl87l |. In brief, this 
method is based on the high probability for annihilation of p on a single peripheral nucleon 
which leads to a residual nucleus containing one neutron or one proton less than the target 
nucleus {Z, N). When the two residual nuclei {Z, N — 1) and {Z — 1, N) are radioactive, then 
measuring their activities can provide the ratio between the probability for p annihilation 
on a neutron to that on a proton. This is based on the reasonable assumption that following 
such extremely peripheral annihilation the resulting pions will not interact with the residual 
nucleus. Considering that the absorption takes place within a narrow range of radii in the 
outer surface region of nuclei, these ratios may provide information on the ratios between 
neutron and proton densities at that region, supplementing the information provided by the 
atomic X-rays. 

Experimental ratios of p absorption on neutrons to absorption on protons were taken 
from Refs. ji^l, [50^. Initial calculations showed that very large contributions to the resulting 

came from ^°^Cd and ^^^Sn and subsequently these two nuclei were excluded from the 
data set. Possible explanations for the problem with these two nuclei in terms of a pp quasi 
bound state are given in Ref. |l88l |. We have therefore used 17 values of absorption ratios 
for the following nuclei: ^^Ca, ^^Ni, ^^Zr, ^oomo, ^^'^^^Ru, ^^^Cd, ^^^Sn, i28,i30rpg^ i44,i54gj^^ 
i48Nd, leoGd, i76Yb, 232Th and ^s^U. 



C. Analyses of antiprotonic atom X-ray data 

Detailed analyses of the results of the PS209 collaboration have been published in a series 
of papers, dedicated each to a particular subset of the data such as neighboring nuclei or 
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FIG. 31: Global best-fit values for zero-range p-nucleus potentials as function of the 
parameter 7 of Eq.([8]) for three shapes of the neutron density 
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FIG. 32: Same as Fig. [3T]but for finite-range potentials with a Gaussian parameter /3=0.9 fm. 
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FIG. 33: Global best-fit values as function of a Gaussian range /3 for the halo shape for p„ with 
7=1.0 fm. 



isotopes of the same element. In several cases it is necessary to take into consideration the 
effects of possible E2 resonances, when energy of a nuclear E2 transition is very close to 
the energy of the atomic transition being studied. In what follows we discuss only global 
fits to the entire data set of 90 points as part of a study of medium- modification of the pN 
interaction [l^, 98 . 

Figure [31] shows the values for the best-fit potential of the type Eq. fl59|) obtained with 
only two adjustable parameters, the real and imaginary parts of Bq. The halo shape for p„ 
yields the lowest value of but the minimum at 7 ^0.5 fm is unacceptable as representing 
the average dependence of r„ — Vp on the neutron excess, as discussed in Sec. IIII D 2[ In 
Fig. [321 are shown similar results for a finite-range version of the potential, obtained with 
Gaussian folding, as given by Eq. f[TU]) . using a range parameter of (3=0.9 fm. The lowest 
is significanty lower than the corresponding value in Fig. [JT]and is obtained for 7 ^0.9 fm, 
which is a most acceptable value, see Sec. IIII and Ref. |l6|]. The FR parameter /3=0.9 fm 
is chosen because over a range of values of 7 a minimum of is obtained for this value of 
(3=0.9 fm, as seen in Fig. [331 This minimum means a pei' point of 2.2 which is quite 
good. The parameters of the potential are Re6o=l-l±0.1 fm, Im6o=l-8±0.1 fm for 7=1.0 fm 
and 6 = —0.035 fm, see Eq. ([8]). These parameters are not qualitatively distinct from the 



parameters obtained recently by Wycech et al. [189| using somewhat different values for (3 
and 7, and also including a p-wave absorptive term in the p optical potential. 

Figure [311 shows results when the isovector parameter bi is varied in the fit in addition to 
the isoscalar parameter 60 • It is found that Re6i is always consistent with zero (not shown) 
whereas the other three parameters vary monotonically with the neutron radius parameter 
7. It is seen that the minimum of is obtained for the same value of 7 ^1.0 fm as before 
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and the quality of fit is not improved. Moreover, at the best-fit point Im6i is consistent with 
zero. Note that a non-zero value for this parameter will be obtained if a very different value 
of 7 is used to represent neutron densities. 

It is interesting to make a few comparisons between the values of the differences be- 
tween rms radii of neutron and of proton distributions implied by the global best-fit value 
of 7=1.0 fm, and differences obtained in de tailed analyses of a small group of p atoms. For 
example, for ^^'^Sn it is found in Ref. [l85| that r„ — rp = 0.08io;o4 fm whereas our global 



expression yields 0.13±0.02 fm for this difference if we assign from Fig. [32] an estimated un- 
certainty of ±0 .1 fm to the slope parameter 7. Likewise for ^^^Sn the values are 0.14±0.03 fm 
from Ref. |l85 | a nd .16=b0.02 fm from the present global analysis. Taking ^osp]^ another 
example, in Ref. [193 the rms difference is 0.16±0.04 fm whereas the present global expres- 
sion leads to 0.18±0.02 fm. Similar agreements are found in other cases. 



D. Radial sensitivity of X-ray data 

Before proceeding to the radiochemical results which provide information on the annihi- 
lation of atomic p at the extreme periphery of the nucleus, it is instructive to examine the 
radial sensitivity of X-ray data in order to get some idea on what are the nuclear regions 
that determine the potentials derived above. Following the preliminary results of the 'notch 
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FIG. 35: Functional derivatives of the global best-fit with respect to relative changes in the full 
complex potential (solid curve), in the real part (dotted) and in the imaginary part (dashed) of 
the potential as function of the radial position rj, see Sec. IIIID 4[ 

test' of Ref. [ij where it was shown that p X-ray data are sensitive to the potential at radii 
well outside of the nuclear surface, we apply here the functional derivative method, as dis- 
cussed in Sec. IIIID 4[ Figure [35] shows the FDs for the best-fit potential with 7=1.0 fm, 
5 = —0.035 fm, a Gaussian range of 0.9 fm and with 6o=l-l+il-8 fm. The first conclusion 
from this figure is the dominance of the imaginary part of the potential as the FD with 
respect to it follows closely the FD with respect to the full complex potential. The other 
clear feature are the radial regions where the bulk of |FD| is found, indicating the regions 
to which the data are sensitive. Strictly speaking, the FD refers to the optical potential and 
owing to the finite-range folding the relevant density regions are shifted to approximately 
0.5 fm smaller radii, well outside of the half-density radius (at ?7=0), peaking between ri=2 
and ri=6 where the densities are well below 10% of the central nuclear density. 



E. Analysis of X-ray and radiochemical data 



In the analysis of radiochemical data we adopt the approach of Refs. 18, 49, 13], namely 
that the method is sensitive to the neutron to proton density ratio close to 2.5 fm outside of 

In terms of the global parameter t] defined 



the half-density radius of the charge density [1911 
above, that corresponds to ?7 ~ 5, which is within the region of sensitivity of the X-ray data 
but slightly shifted towards larger radii, as seen from Fig. |35l It is therefore interesting to 



see if analyses of the radiochemical data lead to results consistent with what is obtained 
from the X-ray data. The experimental ratios of absorption on neutrons to absorption on 
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FIG. 36: Fits to the combined X-ray and radiochemical data (X+C, dashed curves) compared to 
fits to the X-ray data (X, sohd curves) only. 



protons were therefore compared to 



Im(&o &i)/n 
Im(6o - fei)/p 



(60) 



where are the volume integrals of the neutron and proton densities, respectively, over an 
appropriate range. To check sensitivity to the chosen range of integration we have carried 
out the integration either between 2.0 and 3.0 fm or between 2.5 and 3.5 fm outside of 
the half-density radius of the charge density. For the finite-range potential used here the 
folded densities were used. Atomic wavefunctions were not included in the integrals because 
their effect largely cancels out in the ratios. Moreover, we note that choosing the range of 
integration was guided by the conclusions of Ref. (l9l[ | which were based on properties of 
the atomic wavefunctions. With the potential parameter h\ consistent with zero the ratios 
Eq. (!60!) become independent of the parameters of the potential, but they are found to be 
sensitive to values of r„ — or to the parameter 7. Examining the for the radiochemical 



data as function of 7, it was found (19| that the minimum occurred for 7 ^ 1.0 fm, as was the 
case for the X-ray data, if the integration range was 2.5 to 3.5 fm (but not 2 to 3 fm) outside 
of the half-density radius of the charge density. Thi s result confirms in a phenomenological 
way the theoretical conclusion of Wycech et al. jl91 | that most of the absorption takes place 
close to 2.5 fm beyond the charge radius. Note that due to the exponential decrease of the 
densities at such large radii the integrals are dominated by the densities close to the lower 
limit of the range of integration. 
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Combining the results of the radiochemical technique with the X-ray data, fits were made 
to the two kinds of data put together, a total of 107 points. From Fig. [36] showing results of 
fits to this combined data set, in comparison with results from fits to the X-ray only data, it 
is seen that the overall picture is the same in both cases, with larger values of per point 
for the combined data. In particular, with the minimum of for the same value of 7, the 
same conclusions are reached regarding neutron densities. 

Before closing this section we look into the broader perspective of neutron densities in 
nuclei obtained from antiprotonic atoms. The following conclusions may be made from the 
global analyses presented above: 

• The potential parameters depend mostly on r„— rp and not on the shape of the neutron 
densities, although the values do depend on the shape of pn- 

• The favored shape of pn is of the 'halo' type, 

• The rms radii of p„ are given on the average by Eq. ([H]) with the parameter 7 1.0 fm. 

A possible difficulty regarding nuclear densities is that the sensitivity of p atom data is to 
extremely small densities, of the order of 5% of the central nuclear density, where the proton 
densities too are not determined well by the traditional methods of electron scattering and 
muonic X-rays. In particular, the 2pF parameterization need not be appropriate to describe 
the outer reaches of the proton densities pp. On the other hand, the present analyses of p 
atomic data lead only to conclusions on differences between neutrons and protons in N ^ Z 
nuclei, both on the differences of rms radii and on differences in shapes. The preferred 
'halo' shape in this context means that the diffuseness parameter a„ is larger than the 
corresponding parameters for protons, which is quite reasonable considering the binding 
energies of least bound nucleons and effects of the Coulomb potential. The disagreement 
with pionic atoms regarding the shapes of p„ is presumably due to the extreme simplification 
introduced in assuming 2pF parameterizations for the densities. Recall that pionic atom data 
are sensitive to densities up to the full nuclear density. It is, therefore, no wonder that there 
are some differences in conclusions obtained from experiments that are sensitive to different 
density regions of the nucleus. The fact that the rms radii obtained with the two methods 
are in full agreement with each other is not a coincidence. It was shown in Sec. IIII D 21 that 
potential parameters for pions depend mostly on r„ — rp and not on the shape of the neutron 
densities. The same results are found for p atoms, as emphasized above. 

F. Deeply bound antiprotonic atom states 

Much the same as with kaonic atoms, see Sec. IIVE| the optical potential for antiprotonic 
atoms is dominated by its imaginary part which is of the order of 100 MeV deep when 
extrapolated into the full nuclear density. Such an absorption inevitably produces strong 
suppression of the atomic wavefunction inside the nucleus, which for a normalized atomic 
wavefunction located mostly outside of the nucleus is equivalent to repulsion. The role of 
the phenomenological attractive real part of the potential is more difficult to visualize. For 
sufficiently attractive potentials there is the possibility of accommodating strongly bound 
nuclear states, albeit very broad, as their wavefunctions are confined to the nuclear volume. 
Due to the orthogonality requirement of nuclear and atomic wavefunctions having the same 
Z-values, the latter might be shifted considerably by the real potential, as demonstrated 
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FIG. 37: Calculated energies of p atomic states in ^"Zr. The lowest energy for each I value 
corresponds to n = Z + 1. The bars represent the widths of the states. 



for kaonic atoms in Fig. [25l where large repulsion of the atomic wavefunction is observed 
as a result of an attractive strong interaction. It is, therefore, n ot surpris ing, that the 



133. 134 in numerical 



phenomenon of saturation of widths of atomic states is observed 
calculations of antiprotonic atom spectra. 

Figure [37] sh ows a calculated energy spectrum for p atoms of Zr as a typical medium- 
weight nucleus jl34| - The strong-interaction potential is taken from global fits to p atomic 
data. The saturation of the widths is easily seen, with the widths increasing very little when 
/ goes down towards 1=0. 



G. Antiproton-nucleus interaction across threshold 

With rather well-established phenomenology of the interaction of antiprotons with nuclei 
in the subthreshold atomic regime, it is of interest to see if the same picture prevails also 
above t hresh old. Indeed early analyses of elastic scattering of 47 MeV antiprotons on carbon 
showed I92I ] that very good fits to the scattering data and to the then available p atom data 



could be obtained with a common optical potential, dominated by its imaginary part and 
based on Gaussian folding with a range parameter of (3=1.2 fm, quite similar to the present 
results. 

Antiprotons offer a unique tool, compared to other exotic atoms, for studying the in- 
teraction with nuclei very close to, but above threshold, in the form of p annihilation. At 
very low energies, below the pp nn charge-exchange threshold, the total p reaction cross 
section consists only of p annihilation. Therefore measurements of annihilation cross sec- 
tions at such low energies may be compared with total reaction cross sections calculated 
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with the optical potentials obtained from fits to antiprotonic atom data. A measurement 
of antiproton annihilation cross section at 57 M eV/c (1.7 MeV kinetic energy) on Ne was 



reported by Bianconi et al. [451] and was shown [193l | to agree with predictions made with 
potentials obtained from global fits to p atom data. Comparisons for similar measurements 
of p annihilation on ^He showed that predictions of annihilation cross sections made with 
p potentials obtained from fits to atomic data for medium-weight and heavy nuclei, do not 
agree with experiment. In contrast, when parameters of the potential were obtained from 
fits only to p atoms of ^'^He, then the predicted annihilation cross section on ^He agreed with 
the measured one. From this example it may be concluded that the p-nucleus potentials 
cross smoothly the threshold from atomic states to the scattering regime. On the other 
hand the global p potentials which reproduce very well p atomic data for targets heavier 
than A ^10, fail to describe similar data for the He isotopes. 

A special case in this context is the annihilation of p on the proton very close to threshold. 
The pp total annihilation cross section was measured at four momenta between 38 and 
70 MeV/c [i^ and with the availability of strong interaction shift, width and yield for the 
Is and 2p levels in antiprotonic hydrogen it is possible to study the pp interaction across 



threshold [194J |. It is found that a Gaussian potential with a range parameter between 1 and 
2 fm produces very good fits separately to the annihilation cross sections and to the atomic 
pB. data, but if both types of data are to be fitted together then the range parameter turns 
out to be /?=1.5±0.15 fm with Re6o = -0.15±0.15 fm, Im6o=1.80±0.06 fm. It is, therefore, 
possible to cross smoothly the borderline of E=0 also for the pB. system. However, the 
interaction parameters are different from those valid for ^'^He and from those valid for target 
nuclei heavier than A ^10. This demonstrates the limitations of using optical potentials 
down to the very light nuclear targets, where the energy dependence of those NN partial- 
wave amplitudes which may have quasibound states or resonances near threshold needs to 



be considered explicitly |l95 . 196l 



Finally, it is interesting to note that the saturation of widths predicted for antiprotonic 
atom states is also predicted and observed above threshold in the form of saturation of 
reaction cross sections jl94| . There is an interesting analogy between widths of bound states 
and total reaction cross sections where for the Schrodinger equation the latter is given by 

o-i? = -^/lmKpt(r)|^(r)|2rfr , (61) 

where ^^(r) is the p-nucleus elastic scattering wavefunction and v is the cm. velocity. Recall 
that the width of a bound state, as discussed in Sees. IIII CI and IIVEI is given by 

_ JImKpt(r)|^(r)Pdr 

where V'(r) is the p full atomic wavefunction. The modification of this expression for the KG 
equation is mentioned in Sec. lIVEi It is therefore to be expected that large local variations 
of the wavefunction, in both cases, is a common mechanism behind departures from linear 
dependence on the imaginary potential. 

At very low antiproton energies where Coulomb focusing is effective, the annihilation 
cross sections on nuclei are expected to scale as ZA^^^ in the perturbative regime |l94 



but the experimental annihilation cross sections on Ne and ^He diffe r strongly from 



this scaling law. This difference is a manifestation of saturation |l93l . 194], confirming the 
general property of saturation of widths as discussed above. 
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VII. THE REPULSIVE S NUCLEAR POTENTIAL 



A. Preview 



One Boson Exchange (OBE) models fitted to the scarce low-energy YN scattering data 
produce within a G-matrix approach, with one exception (Nijmegen Model F), as r nuch 
attraction for the S nuclear potential as the y do for the A nuclear potential, see Ref. 
for a review of 'old' models and Ref. 1198 



197 



^ ^ for the latest state of the art for Njjmegen 

models. Indeed, the best-fit tofrP potential for E~ atoms was found by Batty et al. |199l . 
to be attractive and absorptive, with central depths for the real and imaginary parts of 
25-30 MeV and 10-15 MeV, respectively. It took almost a full decade, searching for E 
hypernuclear bound states at CERN, KEK and BNL, before it was realized that except 
for a special case for f^He, the observed continuum S hypernuclear spectra indicate a very 
shallow, or even repulsive S nuclear potential, as reviewed by Dover et al. j20l| . These 
indica tion s have received firm support with t he m easurement of several {K~ , tt^) spectra at 

Recently, with measurement s of the E~ 

it has 



BNL lioj followed by calculations for ^Be (203 



spectrum in the (vr , K^) reaction taken at KEK across the periodic table j5l|, l204 



become established that the S nuclear interaction is strongly repulsive. In parallel, analyses 
of S~-atom in the early 1990s, allowing for density dependence or departure from the tp 
prescription, motivated mostly by the precise data for W and Pb led to the conclusion 



that the nuclear interaction of Ss is dominated by repulsion ^206i . i207l . l208j, as reviewed in 



Ref. [ij . This might ha ve in teresting repercussions for the balance of strangeness in the inner 
crust of neutron stars |209l |. primarily by delaying the appearance of S~ hyperons to higher 
densities, if at a ll. The inability of the Nijmegen OBE models, augmented by G-matrix 
calculations (l98l] , to produce S nuclear repulsion is a serious drawback for these mo dels at 



present. This problem apparently persists also in the Juelich model approach |210 |. The 



only theoretical works that provide exception are SU(6) quark-model RGM calculations 



by the Kyoto-Nijata group |211[ |. in which a strong Pauli repulsion appears in th e / = 
3/2, ^Si — ^Di EiV channel, and Kaiser's SU(3) chiral perturbation calculation [212'] which 
yields repulsion of order 60 MeV. 

Below we briefly review and update the S~ atom fits and the recent {n~,K~^) KEK 
results and their analysis. 



B. Density dependent S nuclear potentials from fits to S atoms 



Batty et al. 206, 207l | analyzed the full data set of S atoms, consisting of strong- 



interaction level shifts, widths and yields, in order to constrain the density dependence of 
Vsir). By introducing a phenomenological density dependent (DD) potential of the isoscalar 
form 



Vi](r)~[6o + 5o(p(r)MO))"]p(r) , «>0 , (63) 

and fitting the parameters bo, Bq and a to the data, greatly improved fits to the data are 
obtained. Isovector components are readily included in Eq. fl63l) but are found to have 
a marginal effect. Note, however, that the absorption was assumed to take place only on 
protons. The complex parameter may be identified with the spin-averaged Tj~N scattering 
length. For the best-fit isoscalar potentials, ReVs is attractive at low densities outside the 
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FIG. 38: ReV^pt for DD (solid) and for the geometrical model F (dashed) S~ nuclear potentials 
fitted to S~ atomic data. Vertical bars indicate the half-density radius of the nuclear charge 
distribution. 



nucleus, changing into repulsion in the nuclear surface region. The precise magnitude and 
shape of the repulsive component within the nucleus is not determined by the atomic data. 
The resulting potentials are shown in Fig. [38] (DD, solid lines), where it is worth noting 
that the transition from attraction to repulsion occurs well outside of the nuclear radius, 
hence the occurrence of this transition should be largely model independent. To check this 
last point we have repeated the fits to the atomic data with the 'geometrical model' F of 
Sec. IIV| using separate tp expressions in an internal and an external region, see Eq. (1481) . 
The neutron densities used in the fits were of the skin type, with the r„ — parameter 
Eq. ([8]) 7=1.0 fm. The fits deteriorate significantly if the halo type is used for the neutron 
density. The fit to the data is equally good with this model as with the DD model, (x^ 
per degree of freedom of 0.9 here compared to 1.0 for the DD model) and the potentials 
are shown as the dashed lines in Fig. [381 The half-density radius of the charge distribution 
is indicated in the figure. It is clear that both models show weak attraction at large radii, 
turning into repulsion approximately one fm outside of that radius. 

Further insight into the geometry of the S-nucleus interaction is gained by inspecting 
the functional derivatives (FD) of with respect to the optical potentials, see Sec. IIIID4[ 
Figure [3ni shows the FDs based on the best fit of the geometrical model F as discussed above. 
From the differences between the FD with respect to the full complex potential and the FD 
with respect to the real potential it is concluded that both real and imaginary parts play 
similar roles in the E-nucleus interaction. The bulk of |FD| is in the range of 0.5 < ?7 < 6, 
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FIG. 39: Functional derivatives of with respect to the real (solid) and with respect to the full 
complex (dashed) optical potentials for the best fit F potential. 



covering the radial region where the weak attraction turns into repulsion. Obviously no 
information is obtained from S~ atoms on the interaction inside the nucleus. It is also 
interesting to note quite generally that such potentials do not produc e bo und states, and 
this conclusion is in agreement with the experimental results from BNL j202l | for the absence 
of S hypernuclear peaks beyond He. 

Some semi-theoretic al sup port for this finding of inner repulsion is given by RMF calcu- 
lations by Mares et al. (208l ] who generated the S-nucleus interaction potential in terms of 
scalar (a) and vector {u, p) meson mean field contributions, fitting its coupling constants 
to the relatively accurate S~ atom shift and width data in Si and in Pb. The obtained 
potential fits very well the whole body of data on S~ atoms. This potential, which is 
generally attractive far outside the nucleus, becomes repulsive at the nuclear surface and 
remains so inward in most of the acceptable fits, of order 10-20 MeV. The Pb data j205 
are particularly important in pinning down the isovector component of the potential which 
in this model is sizable and which, for E~, acts against nuclear binding in core nuclei with 
N — Z > 0, countering the attractive Coulomb interaction. On the other hand, for very 
light nuclear cores and perhaps only for A = 4 hypernuclei, this isovector component (Lane 
term) gene rates binding of S"^ configurations. In summary, the more modern fits to S~ 
atom data 206 . 207 . 208l | and the present fits with the geometrical model support the pres- 
ence of a substantial repulsive component in the S-nucleus potential which exclude s nor mal 
E- nuclear binding, except perhaps in very special cases such as |He j213 , 214 . 215 , 216 
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FIG. 40: Inclusive (vr ji^"*") spectra on Ni, In and Bi, fitted by a S-nucleus WS potential with 
depths Vb = 90 MeV, Wq = -40 MeV [si]. 



C. Evidence from (tt , K+) spectra 

A more straightforward information on the nature of the S-nuclear interaction has been 
provided by recent m easu rements of inclusive (vr", fC"*") spectra on medium to heavy nuclear 



targets at KEK [5l|, l20J]. The inclusive (71",/^"*") spectra on Ni, In and Bi are shown in 
Fig. HO] together with a fit using Woods-Saxon potentials with depths Vq = 90 MeV for the 
(repulsive) real part and = —40 MeV for the imaginary part. These and other spectra 
measured on lighter targets suggest that a strongly repulsive S-nucleus potential is required 
to reproduce the shape of the inclusive spectrum, while the sensitivity to the imaginary 
(absorptive) component is secondary. The favored strength of the repulsive potential in this 
analysis is about 100 MeV, of the same order of magnitude reached by the DD E~ atomic fit 
potential shown in Fig. [38] as it 'enters' the nucleus inward. The general level of agreement 
in the fit shown in Fig. SD] is satisfactory, but there seems to be a systematic effect calling 
for more repulsion, the heavier is the target. We conclude that a strong evidence has been 
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FIG. 41: Comparison between DWIA calculations [219] and the measured ^^Si(7r^, A'^) spec- 
trum [5]| using six S-nucleus potentials, (a)-(c) with inner repulsion, (d)-(f) fully attractive. The 
solid and dashed curves denote the inclusive and A conversion cross sections, respectively. Each 
calculated spectrum was normalized by a fraction Z^. The arrows mark the S~ — ■^''Alg.s. threshold 
at u; = 270.75 MeV. 



finally established for the repulsive nature of the S-nucleus potential. 

More sophisticated theoretical analyses of these KEK (7r~,fC+) spectra 217 . 218 . 219l . 

22Cll | have also concluded that the S-nuclear potential is repulsive within the nuclear volume, 



although they yield a weaker repulsion in the range of 1 0-40 MeV. An example of a recent 
analysis of the Si spectrum is shown in Fig. Hi] from Ref. |219 | where six different S-nucleus 



potentials are tested for their ability within the Distorted Wave Impulse Approximation 
(DWIA) to reproduce the measured ^^Si(7r~, fT^) spectrum 511- This particular DWIA 
version was tested on the well understood ^'^Si(7r"'", K^^ quasi-free A hypernuclear spectrum 
also taken at KEK with incoming pions of th e same momentum piab = 1-2 GeV/c. Pot entia l 



(a) is the DD, type A' potential of Ref. (2o3, (b) is one of the RMF potentials of Ref. |208 
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that with = 1, and (c) is a local-density approximation version of a G matrix constructed 
from the Nijmegen model F. These three potentials are repulsive within the nucleus but differ 
considerably there from each other. Potentials (d)-(f) are all attractive within the nucleus, 
with (f) being of a t^sP form. All of the six potentials are attractive outside the nucleus, 
as required by fits to the 'attractive' S~ atomic level shifts. The figure shows clearly, and 
judging by the associated x^/N values, that fully attractive potentials are ruled out by the 
data and that only the 'repulsive ' E-nucleus potentials reproduce the spectrum very well, 
but without giving preference to any of these po tenti als (a)-(c) over the other ones in this 
group. It was shown by Harada and Hirabayashi |220 |. furthermore, that the (tt", K^) data 



on targets with neutron excess, such as ^°^Bi, also lack the sensitivity to confirm the presence 
of a sizable (repulsive for S~) is ovector component of the E nucleus interaction as found in 
the S--atom fits [206l. [207l.[208j. 



VIII. E HYPERONS 

Experiments on exotic atoms of S~ hyperons have not been reported so far but the 
possibilities of conducting such experiments have been discussed by Batty et al. |22l| in 
some detail. Here we summarize only the main features of that study, within the broader 
context of the present Review. 

A. Preview 



Dedicated experiments with stopped S hyperons had been proposed in Refs. |22l.[223l. 



223 in order to produce some of the lightest AA hypernuclei, \^B.e and j^^B. (if the latter 
is particle stable), and /^B, respectively, by looking for a peak in the outgoing neutron 
spectrum in the two-body reaction 

S- + UZ-l) + n. (64) 



These proposals motivated the AGS experiment E885 |225l | on ^^C, using a diamond target to 
stop the S~ hyperons resulting from the quasi-free peak of the p{K~, i^'+)S~ initial reaction. 
An upper bound of a few percent was established for the production of the /^B hyp ernucleus . 



The experimental evidence for ^^He and ^^H had to await different techniques [2261 . 1227 



although the evidence for the latter species remains controversial. The stopped S reaction 
in deuterium, (S~(i)atom Hn, was used in the AGS exp eriment E813 to search for the 



doubly strange H dibaryon, yielding a negative result |228l |. A similar search by the KEK 
E224 collaboration, stopping S~ on a scintillating fiber active carbon target, also yielded a 
nega tive result j229j . On the positive side, following the discovery of a double-A hypernucleus 



230| in light emulsion nuclei by the KEK stopped S experiment El 76, and its interpretation 



due to xaB 231\ . this experiment gav e evidenc e for several events, each showing a decay into 



a pair of known single A hypernuclei j232l . 233l |. One could then attempt to use these events 



in order to deduce properties of the initial S~ atomic states. However, the typical error of 
100 keV incurred in emulsion work is three orders of magnitude larger than the anticipated 
sensitivity of strong-interaction shifts and widths of S~ atomic levels to the S-nucleus strong 
interaction. This simple argument provides a major justification for pursuing a program of 
measuring S~ X rays, in parallel to more conventional strong-interaction reactions involving 
S hyperons, as discussed in Sec. IVIIIBI 
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Very little is established experiment ally or phenomenologically on the interaction of S 
hyperons with nuclei. Dover and Gal [234 1. analyzing old emulsion data which had been 
interpreted as due to S~ hypernuclei, obtained an attractive S-nucleus interaction with a 
nuclear potential well dept h of Vi"''= 21-24 MeV. This range of values agreed well with 
their t heor etical prediction jl97l | for S in nuclear matter, using model D of the Nijmegen 
group j235l | to describe baryon-baryon interact ions in an SU(3) picture, in contrast with the 



H-nucleus repulsion obtained using model F |236| . Similar prediction s we r e su b sequ ently 



made with more detailed G matrix evaluations by Yamamoto et al. [237 . 238 . 239 who 
argued for a considerable A dependence of such that the well depth for light and 

medium weight nuclei is significantly lower than for heavy nuclei where it approaches the 
value calculated for nuclear matter. It should be noted, however, that the predictions of 
the Nijmegen hard-core models D and F are extremely sensitive to the values assumed for 
the hard-core radius. The confidence in the predictive power of model D in strangeness 
—2 hypernuclear physics was to a large extent due to its success to yield the substantial 
attractive AA interaction which was deemed necessary to reproduce the three known AA 
binding energies in the 1990s. This picture has changed during the last decade for several 
reasons, as follows. 



Inclusive {K , K^ ) spe ctra taken at the KEK-PS and at the BNL-AGS accelerators 
on 12 C, Refs. |241| respectively, when fitted near the S -hypernuclear threshold 



yield more moderate values for the attractive S well depth, Vq'"'' ~ 15 MeV. 



The uniquely identified fj^e hypernucleus 226| implies a considerably weaker AA in- 



teraction than produced by reasonable versions of Model D. In particular, the Nijmegen 
soft-core potentials NSC97 [242] provide a more realistic frame work for reproducing 



the weaker strength of the AA interaction, as discussed in Refs. |243l . |244 



New versions of Nijmegen extended soft-core potentials ESC04 198l . 245 ] predict a 
weak S-nucleus interaction with a dehcate pattern of spin and isospin dependence. 
Similar conclusions are also reached in spin-flavor SU(6) quark models by Fujiwara et 

al. rara. 



Looking ahead at the prospects of further research in this strangeness —2 sector, it is 
safe to argue that if the interaction of S hyperons with nuclei is sufficiently attracti ve t o 
cause binding, as has been repeatedly argued since the original work of Dover and Gal 234| , 
then a rich source of spectroscopic information would become available and the properties 
of the in-medium 'EN interaction could be extracted. Boun d states o f S hypern uclei would 
also be useful as a gateway to form double A hypernuclei |248l . |249| . |250| . [251I . Finally, a 



minimum strength for Vo""* o f abo ut 15 MeV is required to realize the exciting possibility 



of 'strange hadronic matter' |252| . where protons, neutrons. As and Ss are held together 
to form a system which is stable against strong-interaction decay. The study of S-nuclear 
interactions, as part of studying strangeness —2 hadronic and nuclear physics, is high on the 
agenda of two forthcoming major high- intensity hadron facilities. 

• At J-PARC, Japan, the main accelerator ring is a 50-GeV proton synchronton and the 
proton beam, with 30 GeV energy and 9/iA current initially, will produce various high- 
intensity beams of secondary particles. Strangeness —2 physics will be explored with a 
K~ beam at piab = 1-8 GeV/c. An approved day-1 experiment is Spectroscopic study of 
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the E-hypernucleus WBe via the ^'^C{K , K^) reaction (T. Nagae, Spokesperson |253 |). 



The overall energy resolution in the S~ bound-state region is expected to be better 
than 3 MeV at FWHM, using an improved version of the existing SKS spectrometer 
at KEK. Another J-PARC approved experiment, although not prio ritiz ed as 'day-1', 



is Measurement of X rays from S atoms (K. Tanida, Spokesperson [253j), the physics 



considerations and the experimental concerns of which are discussed in Sec. IVIIIBI 



below, following the work of Batty et al. [221 



A major component of the upgraded GSI facility in Darmstadt, Germany, will be the 
High Energy Storage Ring (HESR) for high-intensity, phase-space cooled antiprotons 
between 1.5 and 15 GeV/c. A general purpose detector PANDA (Proton ANtiproton 
at DArmstadt) will be set up at the HESR. PANDA is scheduled to provide access 
to high-resolution spectroscopy of 5* = —2 hype rnuclei and hyper-atoms by producing 
abundantly S~ hyperons via the reactions [254 



p + p^H p + + (65) 

occurring on a nuclear target at puh ~ 3 GeV/c. The trigger for these reactions will 
be based on the detection of high-momentum H antihyperons at small angles or of 
mesons produced by the absorption of antihyperons in the primary target nuclei. The 
S~ hyperons will be slowed down and captured in a secondary nuclear target. One 
expects in this way to reconstruct approximately 300 stopped S~ hyperons per day 



in PANDA. A recent simulation is found in Ref. 1255 



B. H atoms 

Conventional measurements of particle energies to investigate S~ hypernuclei suffer from 
insufficient accuracy for providing detailed quantitative information on the interaction of 
S~ hyperons with nuclei. Complementarily, the usual precision for measuring the energies 
of X-rays from transitions between levels of exotic atoms offers the possibility of obtaining 
further information. The successful observation and reasonably precise measurement of 
strong interaction effects in S~ atoms, which had provided significant clues to the interaction 
of E~ hyperons with nuclei, may serve as a guide in assessing the feasibility of experiments 
on exotic atoms of S~. Recall that the H~ and the S~ hyperons have very similar masses 
and lifetimes, namely 1321.32 vs. 1197.34 MeV, and 0.1642 vs. 0.1482 nsec, for S~ and 
S~ r espectively. Full atomic cascade calculations were performed for S~ and S~ atoms 



221| and confirmed that the processes within these two hadronic atoms are very similar. 
The remaining major differences are in the production reactions. Whereas relatively slow 
S~ hyperons are produced by the p(ii'~, 7r+)S~ stopped K~ reaction, the p{K~ , K^)E~ in- 
flight reaction produces relatively fast E~ hyperons, thus causing non-negligible decay losses 
during the slowing down time of the S~ hyperon. Prior to such an experiment it is necessary 
to optimize the experimental setup, which includes a hydrogen production target, a heavy 
moderator such as Pb or W, the target to be studied and the detectors, both for X-rays and 
for the detection of the outgoing K^, which is essential in order to reduce background. 

When selecting targets for possible experiments on S~ atoms, it must be assumed that 
such experiments will probably not be feasible on more than very few targets, and one must 
therefore ask whether it is at all likely that useful information on the interaction of S~ 
with nuclei will be obtained from the resulting rather limited range of data. It was shown 
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50 52 
atomic number 

FIG. 42: Solid curves: calculated strong interaction widths and upper level relative yields for the 
7i level in medium-heavy H~ atoms as function of the atomic number Z. The dashed curves are 
for bo = —0.25 + i0.04 fm, i.e. a repulsive real potential. 



(256 | that the main features of the interaction of K and S with nuclei, as found from 



analyses of all the available data, may in fact, be obtained by analyzing a small fraction 
of the available hadronic atom data, if the target nuclei are carefully selected. A key point 
here is to have target nuclei over as wide a range of the periodic table as possible. This 
observation suggests that experiments on S~ atoms may provide useful information. 

For estimating strong interaction effects in E~ atoms the tp potential Eq. ([3]) was adopted 
with Re6o = 0.25 fm which yields a potential depth of about 20 MeV inside nuclei, and 
Imfeo = 0.04 fm, yielding for the imaginary potential a depth of about 3 MeV. Whereas the 
real potential may be regarded as 'typical', acco rding to the above discussion, the imaginary 
potential is about twice as large as estimated |237l | in model D. Reducing the imaginary 
potential will only cause the calculated widths of the states to decrease by roughly the same 
proportion, and the relative yields (see below) of transitions to become larger. This will 
not, however, change the last observed atomic level. In choosing criteria for the suitability 
of a transition as a source of information on the S nucleus interaction, one is guided by 
experience with other hadronic atoms [1] and select X-ray transitions (n + 1, / + 1) (n, /) 
between circular atomic states (n = / + 1) with energies greater than 100 keV, where the 
strong interaction shift for the 'last' (n, /) level is at least 0.5 keV and the width less than 
about 10 keV. The 'upper' level relative yield, defined as the ratio of the intensity of the 
(n + 1, / + 1) (n, /) X-ray transition to the summed intensity of all X-ray transitions 
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FIG. 43: Calculated strong interaction widths and upper level relative yields for the 4/ level in H~ 
atoms. 

feeding the (n + 1, Z + 1) state, is also required to be at least 10%. 

Strong interaction shifts and widths of S~ atomic levels have been calculated using the 
above optical potential for a large number of nuclei. As the overlap of atomic wavefunctions 
with nuclei vary smoothly with charge number, it is to be expected that generally shifts, 
widths and yields will vary smoothly along the periodic table. Figure |12] shows calculated 
widths and 'upper' level relative yields for the 7i state in medium-heavy S~ atoms and it 
is seen that a suitable target may be found near Sn or I. The dashed lines in this figure are 
obtained by reversing the sign of the real potential used for calculating the solid curves. It is 
seen that in such a case the range of suitable targets will move to between I and Ba, where 
the strong interaction width and relative yield are more acceptable. The sign of the strong 
interaction shift will be reversed in this case, but it has no experimental consequences. This 
exemplifies a general property of hadronic atoms, which are dominated by the Coulomb 
interaction, namely, that large variations in the strong interaction potential will move the 
proposed targets only a few units of charge along the periodic table. 

Figure US] shows results for the 4/ state of atoms, where it is seen that for a Si target 
the effects could be too small to measure, whereas for Ca the width could be too large 
and the relative yield too small. In this region a CI target may be appropriate, perhaps 
in the form of the liquid CCI4. More detailed results for CI are shown in Fig. SH where 
the sensitivities to assumptions regarding the optical potential are also typical of results 
for other targets. The solid curves connect points obtained within the tp potential for fixed 
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FIG. 44: Calculated strong interaction shifts and widths for the 4/ level in H atoms of CI for 
different optical potentials, see text for detail. 



TABLE VIII: Predictions for likely targets for a atoms experiment. Calculations are based on 
a tp potential with 60 = 0.25 + i0.04 fm. E^; is the transition energy, Y is the upper level relative 
yield. 



target 


F 


CI 


Sn 


I 


Pb 


transition 


4/^3(i 


55^4/ 


8j ^ 7i 


8i ^ li 


10/ 9k 


E, (keV) 


131.29 


223.55 


420.25 


474.71 


558.47 


Y 


0.31 


0.37 


0.76 


0.43 


0.58 


shift (keV) 


1.56 


1.84 


0.67 


2.79 


1.73 


width (keV) 


0.99 


1.14 


0.43 


2.21 


1.26 



values of Re6o, listed above the lines. The four points along each line correspond to values of 
Im6o from 0.05 fm down to 0.02 fm in steps of 0.01 fm. Departures from this tp potential are 
represented by the dashed lines, calculated from phenomenological density dependent (DD) 
real potentials similar to those found from analyses of experimental results for K~ and S~ 
atoms, as discussed in Sees. HVland rvTTt respectively. The imaginary part of the potential is 
of the tp type and the points along the dotted lines correspond to the same values of Imfeg 
as above. The real potentials in these calculations are similar to the real potential for S~ 
atoms, having an attractive pocket about 5-10 MeV deep outside the nuclear surface, with 
a repulsive potential of about 20-30 MeV in the nuclear interior. The results in the figure 
serve only to illustrate the expected range of strong interaction effects. If the actual values 
of shift and width turn out to be within the area covered by the lines, these effects will most 
likely be measurable. 
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Table IVIIII summarizes results for possible targets for S atom experiments |221 |. It 



should be kept in mind that due to the discrete nature of quantum numbers it is not always 
possible to 'fine-tune' one's choice of a target, considering widths and yields, in spite of their 
smooth variation with Z. Furthermore, as mentioned above, the main difficulty is likely to 
be associated with the efficient slowing down of S~ hyperons. 
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